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Abstract 


Evaluated  nuclear  data  for  186Re  identify  the  majority  of  spin-parity  assignments 
as  tentative,  with  approximate  values  for  the  energies  of  several  levels  and 
transitions.  In  particular,  the  absence  of  known  transitions  that  feed  the  Jn  =  8+ 
isomer  186mRe  motivates  their  discovery.  This  isomer,  which  has  a  half  life  of 
2  x  105  years,  has  a  potential  application  in  an  isomer  power  source.  Additionally, 
the  isomer’s  role  in  certain  nucleosynthesis  processes  is  not  well  understood,  so 
measured  cross  sections  for  transitions  that  feed  the  isomer  would  have 
astrophysical  implications.  Using  the  GErmanium  Array  for  Neutron  Induced 
Excitations  (GEANIE)  spectrometer  at  the  Los  Alamos  Neutron  Science  Center 
(LANSCE),  (n,  2ny)  and  (n,  n’j)  reactions  in  a  99.52%  enriched  187Re  target  were 
used  to  obtain  y-ray  spectra  from  186 Re  and  18 'Re,  respectively.  The  experimental 
data  reveal  5  new  transitions  in  186Re  and  4  new  transitions  in  18 'Re.  Similarities 
between  the  level  schemes  of  184Re  and  186Re  suggest  that  one  of  the  newly-observed 
transitions  in  186 Re  feeds  the  isomer  from  a  level  at  414.9  keV.  The  y-ray  energy 
measured  for  this  transition  implies  an  isomer  energy  of  148.2(5)  keV,  which  is  a 
significant  improvement  over  the  adopted  value  of  149(7)  keV. 
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ANALYSIS  OF  THE  NUCLEAR  STRUCTURE  OF  186RE  USING 
NEUTRON-INDUCED  REACTIONS 

I.  Introduction 


1.1  Background 

Natural  rhenium  consists  of  37.4%  185Re  and  62.6%  18'Re.  It  was  the  last  stable 
element  to  be  discovered,  and  it  is  notable  for  the  fact  that  it  has  the  highest  boiling 
point  and  the  third  highest  melting  point  of  all  elements.  These  physical  properties 
make  rhenium  suitable  as  a  coating  for  turbine  blades  used  in  aircraft  jet  engines. 
Due  to  its  demand  in  the  aerospace  industry  and  its  rarity,  samples  of  rhenium 
metal  can  be  very  expensive  [1]. 

The  isotope  186Re  can  be  created  by  neutron  capture  (n,  7)  on  185Re  or  by  (n,  2 n) 
reactions  from  18 'Re,  and  is  unstable  in  its  ground  state.  In  the  ground  state,  186Re 
decays  via  electron  capture  or  f3~  decay  with  a  half  life  of  2.718  days,  but  186Re  also 
has  an  unusually  long-lived  metastablc  state,  or  isomer.  The  isomer  186mRe  has  a 
half-life  of  2.0  x  105  years  and  decays  via  7-ray  cascade  to  the  ground  state  [2], 

1.2  Motivation 

The  motivators  for  research  on  the  nuclear  structure  of  186Re  are  (1)  addressing 
a  long-standing  astrophysical  question  related  to  the  rhenium/osmium 
cosmochronometer,  and  (2)  developing  a  radioisotope  power  source  for  use  on  the 
battlefield. 
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Re-Os  Cosmochronometer. 


The  isotope  187Re  has  a  half-life  in  its  ground  state  of  4.35  x  1010  years  [2],  and 
it  f3~  decays  to  the  stable  isotope  18'0s.  The  long  half-life  of  187Re,  on  the  order  of 
the  age  of  the  universe  (15  ±  2  Gyr)  makes  the  187Re/18'Os  system  suitable  as  a 
chronometer  with  which  one  can  estimate  the  ages  of  astronomical  objects. 
Measurement  of  the  relative  abundances  of  18 'Re  and  18 'Os  in  meteorites  permits 
one  to  date  the  nucleosynthesis  of  rhenium  and  osmium  by  high  neutron  flux  events 
such  as  supernovae.  The  Re-Os  radioactive  decay  system  has  been  proposed  as  a 
mechanism  by  which  one  can  determine  the  age  of  the  galaxy  [3] . 

The  nucleosynthesis  of  18 'Re  and  18 'Os  involves  both  a  rapid  r-process,  and  a 
slow  s-process.  In  the  r-process,  neutron  fluxes  in  the  stellar  environment  are  so 
high  that  unstable  isotopes  have  a  high  probability  of  undergoing  neutron  capture 
before  they  have  the  opportunity  to  f3  decay,  while  in  the  s-process  neutron  capture 
rates  are  low  enough  that  (3  decay  plays  an  appreciable  role.  These  processes,  shown 
in  Figure  1,  result  in  the  production  of  18 'Os  almost  entirely  by  s-process  and  18 'Re 
mostly  by  r-process.  In  this  scheme,  186Re  is  considered  a  branch  point,  in  that  it 
can  decay  with  a  half-life  of  3.72  days  into  1860s  or  undergo  neutron  capture  to 
create  187Re  [4,  5,  6].  Consideration  of  the  role  of  the  186mRe  isomer  (which  is  not 
known  to  f3~  decay  independently)  in  this  process  would  reduce  errors  inherent  in 
the  use  of  18'Re/1870s  as  an  effective  chronometer  [7].  Due  to  its  2.0  x  105  year 
half-life,  186mRe  exists  long  enough  to  dramatically  increase  the  chance  of  neutron 
capture  over  /3~  decay,  thus  increasing  the  rate  of  production  of  187Re.  It  is  possible 
that  the  186mRe  isomer  could  be  created  through  neutron-induced  reactions  in  a 
stellar  environment,  which  motivates  the  discovery  of  transitions  that  feed  the 
isomer  and  experimental  determination  of  their  associated  cross  sections. 
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Figure  1.  Rapid  (r,  gray  arrows)  and  slow  (s,  solid  black  line)  processes  involved  in 
the  production  of  18 'Re  and  18 'Os.  Long-lived  nuclides  are  identified  with  bold  boxes, 
with  the  dashed  lines  indicating  weak,  or  secondary  processes.  The  red  outline  around 
the  isomer  186mRe  and  associated  pathways  highlights  the  fact  that  its  role  in  these 
processes  is  not  yet  understood. 


Isomer  Power  Source. 

Long-lived  nuclear  isomers  store  up  to  106  times  more  energy  per  gram  of 
material  than  chemical  energy  sources.  These  isomers  release  their  excess  energy 
through  a  variety  of  decay  mechanisms,  with  half-lives  ranging  from  milliseconds  to 
thousands  of  years.  If  the  isomer  decays  directly  to  the  ground  state,  the  released 
energy  is  not  a  product  of  the  binding  energy  of  the  nucleus,  so  it  does  not  involve 
transmutation.  As  a  result,  the  significant  amounts  of  radioactive  waste  produced 
by  nuclear  fission  in  the  form  of  fission  fragments  and  neutron  activation  products 
are  not  present  following  the  decay  of  nuclear  isomers.  Several  methods,  including 
X-ray  irradiation  and  neutron  bombardment,  have  been  shown  to  induce  depletion 
of  certain  isomers  as  a  means  of  accessing  their  stored  excitation  energy.  However, 
the  goal  of  inducing  a  nucleus  to  release  energy  on  demand  in  practical  quantities 
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remains  elusive  [8]. 

As  technology  on  the  battlefield  becomes  more  advanced,  forward- deployed 
military  forces  require  increasingly  powerful  energy  sources  to  operate  equipment. 
The  portability  of  these  power  sources  can  have  a  critical  effect  on  mission 
readiness,  so  high  energy  density  sources  such  as  nuclear  isomers  are  of  particular 
interest.  To  highlight  the  importance  of  this  goal,  the  Defense  Science  Board  in  its 
October  2013  report  recommended  that  $25  million  be  allocated  to  research  relating 
to  the  development  of  a  radioisotope  power  source  [9] 

i86mjte^  given  its  half-life  of  2.0  x  105  years,  has  the  advantage  of  being  a  stable 
energy  storage  medium.  If  186mRe  could  be  transformed  from  the  isomer  to  the 
ground  state  of  186Re  on  demand,  the  ground  state  would  then  f3~  decay  to  1860s 
with  a  half-life  of  3.7  days,  producing  electrons  with  energies  up  to  1.08  MeV.  In 
one  possible  radioisotope  power  source  design,  the  energetic  electrons  could  produce 
scintillation  light  in  a  crystal  that  could  in  turn  generate  an  electric  current  using  a 
photovoltaic  cell. 

Identifying  transitions  in  the  nuclear  structure  of  186Re  that  feed  the  isomer 
could  lead  to  the  discovery  of  a  method  for  inducing  depletion  of  186,nRe.  The 
potential  of  nuclear  isomers  as  a  power  source  was  outlined  by  P.  Walker  and  J.  J. 
Carroll,  in  which  photon-induced  depletion  of  isomers  was  proposed  as  a  means  of 
achieving  on-demand  energy  release  [8].  The  ultimate  goal  of  these  types  of  isomer 
depletion  experiments  is  the  discovery  of  a  mechanism  by  which  the  isomeric  state 
can  be  depleted  with  a  favorable  energy  balance. 

1.3  Problem 

The  objective  of  this  research  was  to  conduct  an  examination  of  the  nuclear  level 
structure  of  186Re,  with  the  specific  goal  of  discovering  transitions  that  feed  the 
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149(7)  keV  isomer.  186Re  is  an  isotope  with  a  relatively  poorly  known  structure.  A 
search  of  the  Evaluated  Nuclear  Structure  Data  File  (ENSDF)  for  186Re  identifies 
the  majority  of  spin-parity  assignments  as  tentative,  with  approximate  values 
associated  with  the  energies  of  several  levels  and  transitions,  as  shown  in  Figure  2. 
Of  note  in  Figure  2  is  the  the  absence  of  known  transitions  that  populate  the  186mRe 
isomer. 


<0.2  ns 
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0.20  ns 
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Figure  2.  An  extract  from  the  low-energy  level  scheme  of  186Re  shows  that  there  are 
numerous  levels  for  which  the  multipolarity  is  identified  with  parentheses  as  tentative, 
while  transitions  and  energies  assigned  to  several  levels  are  approximate  [10]. 


1.4  Hypothesis 

Measurements  of  y-ray  spectra  from  the  reaction  187Re(n,  2ny)186Re  would 
provide  valuable  input  to  nuclear  structure  databases  such  as  the  Experimental 
Unevaluated  Nuclear  Data  List  (XUNDL)  and  the  Evaluated  Nuclear  Structure 
Data  File  (ENSDF).  New  levels  and  transitions  discovered  in  186Re,  and  in 
particular  levels  that  directly  populate  the  186mRe  isomer,  would  help  solve  the 
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problems  identified  in  Section  1.2. 


1.5  Methods 

Using  the  GErmanium  Array  for  Neutron  Induced  Excitations  (GEANIE) 
detector  array  at  the  Los  Alamos  Neutron  Science  Center  (LANSCE)  Weapons 
Neutron  Research  (WNR)  facility,  7-ray  spectra  obtained  from  neutron- induced 
reactions  were  used  to  identify  new  7-ray  transitions  in  186Re  and  verify  transitions 
and  levels  already  described  in  existing  literature.  The  primary  neutron-induced 
reaction  studied  was  18'Re(n,  2n7)186Re,  from  irradiating  a  rhenium  powder  sample 
enriched  in  18 'Re.  The  isotope  184Re  has  a  similar  structure  to  186 Re,  so  known 
7-ray  transitions  in  184Re  were  used  as  a  guide  to  identifying  transitions  feeding  the 
isomer  in  186Re. 


6 


II.  Theory 


The  physics  of  nuclear  structure,  which  is  concerned  with  the  levels  and 
transitions  in  nuclei  at  excitation  energies  below  approximately  2-3  MeV, 
incorporates  both  collective  and  microscopic  models  of  nucleon  motion  [11],  Both 
are  important  to  understanding  the  structure  of  nuclei  with  odd  numbers  of  protons 
and  neutrons,  such  as  186 Re. 

2.1  Nuclear  Landscape 

The  observed  instability  of  radioisotopes  provided  the  first  clues  to 
understanding  the  physics  underlying  nuclear  structure.  In  a  plot  of  the  atomic 
number  Z  versus  the  neutron  number  N  of  known  isotopes  (Figure  3),  a  clear  valley 
of  stable  nuclei  emerges,  with  unstable  neutron-rich  isotopes  below  the  valley  and 
unstable  proton-rich  isotopes  above  it.  The  shape  of  this  valley  led  to  the 
development  of  the  first  successful  macroscopic,  or  collective,  model  of  nuclear 
structure,  known  as  the  liquid  drop  model  (discussed  further  in  Section  2.3). 

Certain  nuclei,  with  values  of  N  or  Z  equal  to  the  so-called  magic  numbers,  are 
known  to  be  especially  stable,  a  behavior  not  predicted  by  collective  models.  The 
magic  numbers,  N,  Z  —  2,  8,  20,  28, 50,  82, 126,  are  shown,  except  for  N,  Z  —  2,  in 
Figure  3.  Gaps  are  also  observed  at  N,  Z  =  40  (not  identified  in  the  figure)  and 
Z  =  64  for  certain  N  [13].  Evidence  of  the  stability  of  nuclei  with  magic  N  or  Z 
values  led  to  the  development  of  the  nuclear  shell  model,  a  microscopic  model  that 
treats  the  neutrons  and  protons  in  the  nucleus  individually  instead  of  collectively. 

The  shell  model  assumes  a  spherical  nuclear  potential,  limiting  its  applicability 
in  some  cases  where  the  nucleus  has  a  non-spherical  shape.  This  modification  to  the 
shell  model,  known  as  the  Nilsson  model,  takes  into  account  the  deformation  of  the 
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Figure  3.  The  chart  of  the  nuclides,  showing  the  rough  valley  of  stable  nuclei  (black 
dots)  in  the  center  of  the  plot  of  proton  number  Z  versus  neutron  number  TV.  Magic 
numbers  are  identified  by  double  lines,  while  realms  in  which  nuclear  deformation  is 
significant  are  outlined  in  purple  arcs.  The  most  common  locations  in  which  spin-trap 
isomers  are  found  are  highlighted  in  green,  while  K- trap  isomers  are  generally  found 
in  the  area  highlighted  in  orange.  Figure  copyright  1994,  IOP  Publishing.  Reproduced 
with  permission  from  [12]. 


nucleus  in  predicting  nuclear  behavior.  The  realms  in  which  nuclear  deformation 
effects  are  significant  are  shown  as  circles  or  arcs  in  Figure  3.  Of  note,  186Re,  the 
nucleus  studied  in  this  work,  falls  well  within  the  mid-shell  region  with 
82  <  TV  <  126  and  50  <  Z  <  82,  where  deformation  is  a  contributing  factor. 


2.2  Microscopic  Models 


Microscopic  models  of  the  nucleus  are  concerned  with  the  behavior  of  individual 
nucleons  and  their  effects  on  the  stability  of  the  nucleus.  A  number  of  microscopic 
models  exist,  the  most  successful  of  which  has  been  the  nuclear  shell  model  [13]. 


Shell  Model. 


The  nuclear  shell  model  (with  spin-orbit  coupling)  provides  a  basis  for 
understanding  the  behavior  of  nuclear  excited  states.  In  this  basic  model,  which 
assumes  a  spherical  nuclear  potential,  protons  and  neutrons  are  assumed  to  be 
independent  particles,  each  occupying  distinct  shells  in  accordance  with  the  Pauli 
exclusion  principle.  When  spin-orbit  coupling  is  included  in  the  model,  the  resulting 
predictions  of  nuclear  stability  are  in  relatively  close  agreement  with  experiment.  In 
particular,  the  shell  model  exactly  predicts  the  experimentally  observed  magic 
numbers,  and  it  also  predicts  that  N,  Z  =  184  is  magic  in  the  absence  of  empirical 
evidence  [14].  The  individual  particle  levels  predicted  by  the  shell  model  are  shown 
in  Figure  4,  in  which  each  level  is  identified  by  its  angular  momentum  quantum 
number  £,  with  s,p,  d, . . .  referring  to  £  =  0, 1,  2, . . ..  The  capacity  of  a  particular 
level  is  shown  to  the  right  of  the  level,  and  to  the  far  right  is  the  cumulative  number 
of  nucleons  in  that  level  and  all  lower  levels. 

The  shell  model  can  be  used  to  explain  behavior  seen  in  nuclear  level  schemes. 
As  an  example,  consider  nuclei  with  both  N  and  Z  magic.  Level  schemes  for  these 
doubly  magic  nuclei  show  a  large  energy  gap  between  the  ground  state  and  the  first 
excited  state,  which  is  explained  by  the  fact  that  a  significant  amount  of  energy  is 
required  to  promote  a  nucleon  out  of  a  closed  shell  to  the  next  (empty)  shell. 
Application  of  the  shell  model  to  the  case  of  nuclei  away  from  closed  shells  requires 
the  independent  particle  approximation,  discussed  in  greater  detail  in  Section  2.4. 

Transition  rates  between  states  reflect  the  stability  of  each  state  and  are 
governed  by  quantum  mechanical  selection  rules.  Gamma  rays,  which  result  from 
transitions  between  energy  levels,  are  electromagnetic  radiation  with  characteristic 
multipolarity  (i.e.,  Ml,  E 2,  etc.).  Quantum  mechanical  calculations  of  the 
transition  probabilities  include  the  following  selection  rules  for  transitions  between 
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Figure  4.  Energy  levels  of  the  nuclear  shell  model  (including  spin-orbit  interaction). 
The  magic  numbers,  corresponding  to  completely  filled  nuclear  shells,  are  identified 
with  circles.  Reproduced  with  permission  from  [14]. 


energy  levels: 


\  <L<£i  +  £f 

An  =  -1  — >  even  electric,  odd  magnetic  (1) 

Att  =  +1  — *  odd  electric,  even  magnetic. 


Here  L  is  the  multipole  order  (L  >  0)  of  the  electromagnetic  transition,  li  and  £/ 
are  the  angular  momentum  quantum  numbers  of  the  initial  and  final  states.  The 
change  in  parity  An  between  initial  and  final  states  is  equal  to  +1  if  there  is  no 
change  in  parity  and  —1  for  a  change  in  parity.  For  given  initial  and  final  states, 
multiple  transition  multipolarities  may  be  possible,  in  which  case  the  most  likely 
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transition  is  the  one  with  the  lowest  multipole  order.  For  a  given  multipole  order, 
electric  transitions  are  more  likely  than  magnetic  ones.  These  effects  are  revealed 
quantitatively  by  the  Weisskopf  estimates. 

Weisskopf  Estimates. 

The  Weisskopf  estimates,  which  arise  from  calculations  of  radial  transition 
matrix  elements,  provide  a  first-order  approximation  of  the  transition  rates  between 
nuclear  states  [11,  14].  These  estimates  are  based  on  the  assumption  of  a 
particularly  simple  spherical  potential,  so  they  are  directly  applicable  only  to  the 
case  of  spherical  nuclear  states.  Even  in  this  limited  case  they  predict  rates  that  can 
result  in  order  of  magnitude  differences  from  measured  decay  rates. 

Given  the  atomic  number  A  and  the  observed  transition  energy  E  in  MeV,  the 
Weisskopf  estimates  give  the  transition  rate  A  in  units  of  s_1.  For  electric  transitions 
El  through  EA  (dipole  through  hexadecapole),  they  are: 

X(E1)  =  1.0  x  10 14A2/3E3 

X(E2)  =  7.3  x  10 7A4/3E5  (2) 

X(E3)  =  3AA2E7 
X(EA)  =  1.1  x  10"5A18/3F;9. 

For  magnetic  transitions,  the  Weisskopf  estimates  for  the  M 1  through  M 4 
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transition  rates  are: 


A  (Ml)  =  5.6  x  1013£3 

A  (M2)  =  3.5  x  10  7A2/3E5  (3) 

A(M3)  =  16  Aa/3E7 
A(M4)  =  4.5  x  10“6A2£9. 

It  is  important  to  note  that  the  leading  constant  in  these  estimates  decreases  by 
approximately  6  orders  of  magnitude  for  each  incremental  increase  in  multipole 
order.  The  transition  rate  also  depends  on  the  energy  of  the  transition.  As 
transition  energy  E  increases,  the  transition  rate  increases  as  E2L+1,  where  L  is  the 
multipole  order  of  the  transition  [11,  14]. 

Deformed  Shell  Model. 

The  shell  model  has  proven  itself  an  incredibly  powerful  tool  with  which  one  can 
explain  many  aspects  of  nuclear  behavior,  including  excited  state  energy  levels  and 
transition  rates.  These  predictions  agree  relatively  well  with  experiment  for  nuclei 
with  N  or  Z  close  to  the  magic  numbers,  where  nuclei  are  spherical  in  their  ground 
state.  However,  as  mentioned  in  Section  2.1,  the  shell  model  applies  directly  only 
when  the  nuclear  potential  can  be  reasonably  approximated  as  spherical.  Nuclei  in 
the  mid-shell  range  with  mass  150  <  A  <  190  are  generally  non-spherical,  so  the 
spherical  potential  approximation  used  in  the  derivation  of  the  nuclear  shell  model 
is  flawed. 

The  two  symmetries  encountered  in  deformed  nuclei  are  oblate  (frisbee)  and 
prolate  (football)  shapes,  with  the  up  direction  aligned  with  the  axis  of  rotation  (by 
convention,  the  z-axis  in  the  center-of-mass  reference  frame).  That  is,  an  oblate 
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nucleus  extends  further  in  the  equatorial  direction  than  in  the  polar  direction,  and 
the  opposite  is  true  for  a  prolate  nucleus. 

For  a  purely  spherical  nucleus,  the  nuclear  radius  R  is  given  by 

R  =  R0A1/3,  (4) 


where  A  =  N  +  Z  is  the  mass  number  and  Rq  =  1.2  fm.  For  a  non-spherical  nucleus 
with  a  static  quadrupole  deformation,  the  nuclear  radius  becomes 


R 


RoA1/3 


i  +  ,  4>) 


(5) 


where  the  are  the  (quadrupole)  expansion  coefficients  for  the  spherical 
harmonics  Yi M(0,  0)  [11,  13].  In  the  center-of-mass  reference  frame,  the  expansion 
coefficients  «0,  a ±2  can  be  written  in  terms  of  a  deformation  parameter  f3  and  an 
angle  7  as  ay  =  /^cosy  and  a±2  =  4^/3  sin  7.  The  parameter  /3  specihes  the  shape  of 
the  distortion  of  the  nucleus,  with  <  0  for  an  oblate  nucleus  and  (3  >  0  for  a 
prolate  nucleus.  The  angle  7  is  a  measure  of  the  axial  asymmetry  of  the  nucleus, 
ranging  from  7  =  0°  for  an  axially  symmetric  nucleus  to  7  =  30°,  which  corresponds 
to  the  case  of  maximum  axial  asymmetry  [11,  13].  The  diagram  in  Figure  5  shows 
the  equipotential  surfaces  for  each  of  the  shapes  described  above. 
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Figure  5.  Characteristic  shapes  of  equipotential  surfaces  in  deformed  nuclei.  The 
spherical  shape  has  (3  —  0,  while  the  oblate  and  prolate  shapes  have  (3  <  0  and  (3  >  0, 
respectively.  The  axial  symmetry  is  a  function  of  0°  <  7  <  30°,  and  the  end  views  show 
different  degrees  of  axial  asymmetry  for  the  prolate  and  oblate  shapes.  Reproduced 
with  permission  from  [13]. 


A  more  accurate  model  of  the  nucleus  is  based  on  this  deformed  potential,  in 
which  the  degree  of  deformation  e  is  accounted  for  in  determining  the  energy  levels. 
The  degree  of  deformation  e  is  a  function  of  the  deformation  parameter  /3  according 
to  the  relation 


0  = 


(6) 


The  deformed  shell  model,  also  known  as  the  Nilsson  model,  accounts  for  the 
non-spherical  nature  of  mid-shell  nuclei  [2],  The  resulting  energy  levels  show 
significant  dependence  on  the  deformation,  which  is  evident  in  the  so-called  Nilsson 
diagram  of  Figure  6.  For  186Re,  f3  =  0.22,  which  corresponds  to  e  =  0.194  per 
Equation  6  [15]. 

In  the  Nilsson  model,  the  energy  levels  of  deformed  nuclei  depend  on  the 
projection  of  the  angular  momentum  on  the  axis  of  symmetry,  defined  as  the 
quantum  number  K  as  shown  in  Figure  7.  Levels  in  the  Nilsson  model  are  identified 
according  to  the  convention  K77  [NnzA]  .  The  first  quantum  numbers  Iin  define  the 
K  value  and  parity  of  the  state,  N  represents  the  principle  quantum  number  which 
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(a)  Neutron  energy  levels  (b)  Proton  energy  levels 

Figure  6.  Deformed  shell  model  (Nilsson)  energy  level  diagrams  for  (a)  neutrons  and 
(b)  protons  in  deformed  nuclei  with  82  <  N  <  126  and  50  <  Z  <  82.  The  parameter 
e  (e2  in  the  figure)  is  a  measure  of  the  deformation  of  the  nucleus:  e  <  0  for  oblate 
nuclei  and  e  >  0  for  prolate  nuclei.  The  vertical  red  line  identifies  the  value  of  e  =  0.195 
for  186Re.  Magic  numbers  are  identified  by  circles.  Figure  copyright  2009  Wiley-VCH 
Verlag  GmbH  &  Co.  KGaA.  Reproduced  with  permission  from  [2]. 

defines  the  major  shell,  and  nz  is  the  number  of  modes  of  the  wavefnnction  in  the  z 
direction.  The  quantum  number  A  is  the  component  of  the  orbital  angular 
momentum  along  the  symmetry  axis.  Using  these  definitions,  for  a  nucleon  with 
spin  1/2  and  a  spin  projection  £  =  ±1/2  on  the  nuclear  symmetry  axis,  the  K  value 
is  equal  to  A  ±  £. 
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Figure  7.  Definition  of  quantum  number  K  and  angle  of  inclination  0  for  a  single 
nucleon  of  spin  j  orbiting  a  prolate  deformed  nucleus.  Reproduced  with  permission 
from  [13]. 


In  a  deformed  nucleus,  the  single  particle  orbits  can  have  different  inclinations 
relative  to  the  axis  of  symmetry.  The  inclination  affects  the  energy  of  the  particular 
orbit,  depending  on  the  average  distance  between  the  orbit  and  the  rest  of  the 
nucleus.  When  the  orbiting  nucleon  remains  close  to  the  central  nuclear  mass,  it  will 
be  more  tightly  bound  (due  to  the  attractive  strong  nuclear  force)  and  have  a  lower 
energy  than  a  nucleon  that  orbits  at  a  higher  inclination  angle  and  spends  the 
majority  of  the  orbit  at  a  greater  distance  from  the  central  mass.  These  two  orbits 
are  shown  in  Figure  8. 


Figure  8.  Two  nucleon  orbits  for  a  prolate  nucleus  at  different  inclination  angles  6.  In 
this  figure,  the  orbit  labeled  K1  has  a  lower  energy  than  the  one  labeled  K2.  Repro¬ 
duced  with  permission  from  [13]. 


2.3  Collective  Models 


Collective  models  treat  the  nucleus  as  a  whole,  instead  of  accounting  for  the 
interactions  between  individual  nucleons.  The  average  motion  of  the  protons  and 
neutrons  is  used  to  explain  empirical  observations. 

Liquid  Drop  Model. 

The  most  familiar  collective  model  of  the  nucleus  is  the  so-called  liquid  drop 
model.  The  liquid  drop  model  was  initially  developed  by  H.  Bethe  and  C. 
Weizsacker,  and  it  shows  good  qualitative  agreement  with  experiment  in  the  case  of 
predicting  nuclear  binding  energies. 

Bohr-Mottelson  Model. 

Much  like  in  atomic  and  molecular  physics,  in  which  molecular  rotational  and 
vibrational  states  are  evident  from  features  in  visible  or  infrared  spectra,  the  nucleus 
can  exhibit  collective  vibration  and  rotation.  This  collective  motion  is  the  basis  for 
the  nuclear  model  proposed  by  A.  Bohr  and  B.  Mottelson  in  the  1950s. 

The  collective  model  of  Bohr  and  Mottelson  incorporates  asymmetries  in  the 
nucleus,  and  so  accounts  for  the  non-spherical  nature  of  mid-shell  nuclei.  In  their 
model,  in  addition  to  rotations  perpendicular  to  the  axis  of  symmetry  (which 
approximates  a  quantum  rigid  rotor)  the  nucleus  can  experience  (5  or  7  vibrations 
corresponding  to  the  parameters  (3  and  7  described  in  Section  2.2  [11].  In  a  f3 
vibration  the  nucleus  extends  and  contracts  parallel  to  the  axis  of  symmetry  (the 
football  or  frisbee  gets  longer  or  wider,  respectively).  In  a  7  vibration  the  nucleus 
extends  and  contracts  in  a  direction  perpendicular  to  the  axial  direction.  These 
vibrations  are  illustrated  in  Figure  9. 

The  spectral  features  that  the  model  predicts  are  most  clearly  evident  in  7-ray 
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Figure  9.  (3  and  7  vibrational  modes  of  deformed  nuclei.  A  (3  vibration  extends  the 
nucleus  along  its  principal  axis,  while  a  7  vibration  is  a  periodic  variation  in  the  degree 
of  axial  asymmetry.  Reproduced  with  permission  from  [13]. 

spectra  of  even-even  nuclei.  In  these  spectra,  clear  periodicities  are  evident  in  the 
7-ray  energies,  corresponding  to  vibrational  or  rotational  excitation  of  the  nucleus. 
These  can  be  organized  into  distinct  bands,  where  7-ray  transitions  between  bands 
are  generally  hindered  by  quantum  mechanical  selection  rules. 

Application  of  the  collective  model  to  determine  the  vibrational  and  rotational 
energies  of  odd-odd,  odd-even  or  even-odd  nuclei  requires  additional  approximations 
[11,  13].  For  this  purpose,  combined  models  that  incorporate  both  collective  and 
microscopic  approaches  provide  a  means  of  describing  the  structure  of  nuclei  such  as 
186 Re  from  a  theoretical  standpoint. 

2.4  Combined  Models 

An  explanation  of  the  structure  of  nuclei  with  an  odd  number  of  neutrons  or 
protons  requires  a  combination  of  microscopic  and  collective  models  discussed  in 
Sections  2.2  and  2.3.  Combined  models  are  generally  built  on  a  “core  plus  valence 
nucleon”  description  to  explain  nuclear  behavior  [11], 

Independent  Particle  Approximation. 

The  most  basic  combined  model  is  the  independent  particle  approximation,  in 
which  the  nuclear  properties  are  assumed  to  arise  from  the  motion  of  the  unpaired 
valence  nucleons.  This  is  basically  equivalent  to  the  statement  that  the  individual 
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nucleon-nucleon  interactions  of  the  microscopic  model  can  be  approximated  by  a 
common  nuclear  potential,  i.e.  a  mean-field  approach.  This  approximation 
drastically  simplifies  shell  model  calculations,  since  a  difficult  iV-body  problem 
reduces  to  a  problem  in  perturbation  theory.  This  allows  the  shell  model  to  be 
reasonably  applied  to  nuclei  with  N  or  Z  away  from  the  magic  numbers  [13]. 

An  illustrative  example  of  the  independent  particle  approximation  is  the 
estimation  of  the  ground-state  spin  J  of  a  nucleus  with  an  odd  number  of  protons  or 
neutrons.  Since  the  spins  of  paired  nucleons  cancel  to  0,  the  total  spin  J  of  a 
nucleus  with  an  odd  number  of  protons  or  neutrons  will  equal  the  spin  j  of  its  lone 
unpaired  nucleon.  For  an  odd-odd  nucleus,  the  total  spin  is  determined  by  taking 
the  vector  sum  of  the  spins  of  the  two  unpaired  nucleons  [14]. 

The  independent  particle  approximation  can  help  to  explain  why  level  densities 
observed  at  low  excitation  energies  of  even-odd,  odd-even  or  odd-odd  nuclei  are 
greater  than  those  predicted  by  the  shell  model  of  the  nucleus.  In  even-odd, 
odd-even  or  odd-odd  nuclei,  nucleon  spins  (proton,  neutron,  or  both)  are  inherently 
unpaired.  The  energy  required  to  excite  the  nucleus  to  a  higher-energy  state  is 
greater  when  the  nucleus  consists  of  paired  nuclei,  as  energy  is  required  to  break 
apart  nucleon  pairs.  Since  less  energy  is  required  to  excite  nuclei  with  unpaired 
nucleons,  even-odd,  odd-even  and  odd-odd  nuclei  are  found  to  have  higher  level 
densities  at  low  excitation  energies  [14]. 

2.5  Metastable  Excited  States  (Isomers) 

The  long  half  lives  of  nuclear  isomers  are  a  direct  result  of  selection  rules  that 
govern  7  decay.  The  Weisskopf  estimates  are  helpful  for  explaining  the  behavior  of 
nuclear  isomers.  Long-lived  isomers  typically  have  spin-parity  (Jn)  assignments  that 
differ  significantly  from  that  of  the  states  to  which  they  can  decay.  Per  the  selection 
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rules,  a  direct  transition  from  the  isomer  to  a  lower-energy  state  would  have  a  high 
multipole  order.  Thus,  in  a  so-called  spin-trap  nuclear  isomer,  significant  differences 
in  angular  momentum  between  the  initial  and  final  states  result  in  very  low 
transition  rates,  or  long  half  lives  [8,  12,  16].  The  Weisskopf  estimates  also  imply 
that  low-energy  transitions  are  inhibited. 

Another  factor  adding  to  the  long  half  lives  of  certain  isomeric  states  is 
A'-hindrance,  which  is  an  effect  only  observed  in  axially  symmetric  deformed  nuclei, 
generally  in  the  mass-180  region  as  shown  in  Figure  3.  The  orientation  of  the 
nuclear  angular  momentum  J  relative  to  the  axis  of  symmetry,  which  is  given  by  the 
quantum  number  K  as  described  in  Section  2.2,  governs  the  ability  of  the  nucleus  to 
decay  from  an  excited  state  to  the  ground  state.  Selection  rules  derived  from  the 
Nilsson  model  require  that  the  multipole  order  L  be  equal  to  or  greater  than  the 
change  in  K  between  the  initial  and  final  (ground)  state.  This  selection  rule  is  not 
strictly  observed,  and  symmetry  breaking  in  the  model  can  lead  to  transitions  in 
which  L  <  A K,  though  such  transitions  are  hindered.  Isomers  in  which  this  is 
found  to  be  the  case  are  known  as  K- trap  isomers  [12,  16]. 

The  nuclear  shapes  of  spin-trap  and  Jl-trap  isomers  are  illustrated  in  Figure  10. 
Per  the  independent  particle  approximation,  the  angular  momentum  arises  from  the 
motion  of  a  small  number  of  valence  nucleons,  shown  orbiting  the  nuclei  at  the  top 
of  each  panel.  In  the  A'-trap  isomer  shown  in  the  figure,  the  nucleus  in  the  isomer 
state  is  shown  with  its  angular  momentum  projection  90°  to  that  of  the  ground 
state,  such  that  the  AK  between  the  isomer  and  ground  state  is  maximized. 

Population  and  Depletion  of  Isomers. 

Because  of  the  low  probabilities  for  transition  directly  between  the  isomer  state 
and  lower-energy  states,  population  or  depletion  of  nuclear  isomers  is  generally 
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Figure  10.  Schematic  plot  of  excitation  energy  versus  angular  momentum  J  in  spin  trap 
isomers  (a)  and  angular  momentum  projection  K  in  K- trap  isomers  (b).  The  variation 
in  the  curve  shown  in  (b)  is  due  to  a  wobbling  motion  of  the  deformed  nucleus.  In  both 
cases,  the  excitation  energy  has  a  local  minimum  at  the  energy  of  the  isomer  state. 
Figure  copyright  1994,  IOP  Publishing.  Reproduced  with  permission  from  [12]. 


achieved  by  exciting  the  nucleus  to  a  higher-lying  intermediate  states  (IS)  from 
which  it  can  then  decay  to  the  isomer  or  lower-energy  state,  respectively.  These 
intermediate  states  are  more  closely  matched  in  spin  to  both  the  isomer  and 
lower- lying  state  [8].  In  the  population  of  isomers  by  photon  absorption,  a  photon 
closely  matched  in  energy  to  the  energy  of  the  IS  can  resonantly  excite  the  nucleus 
into  the  IS.  From  the  IS,  the  nucleus  can  decay  via  y-ray  cascade  or  internal 
conversion  to  feed  the  isomer,  or  it  can  decay  back  to  the  initial  state  [17,  18].  In 
reverse,  depletion  of  nuclear  isomers  involves  excitation  of  the  isomer  to  an  IS  from 
which  it  can  decay  to  a  lower-energy  state  and  release  the  excitation  energy  of  the 
isomer  in  the  process.  These  two  processes  are  illustrated  in  Figure  11. 

If  a  given  IS  at  energy  E  has  a  transition  width  F,  then  a  photon  of  energy 
E  ±  T/2  has  a  high  probability  of  resonantly  exciting  the  nucleus  into  the  IS.  Each 
intermediate  state  that  feeds  the  isomer  thus  has  a  cross  section  for  resonant  photon 
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Figure  11.  Hypothetical  isomer  depletion  (a)  and  photoexcitation  (b)  schemes.  The 
grey  arrows  identify  resonant  photon  absorption  and  excitation  from  the  initial  state 
(isomer  in  (a),  ground  state  in  (b))  to  the  intermediate  state  (IS),  and  the  dashed 
line  is  the  7-ray  cascade  to  the  final  state  (ground  state  in  (a),  isomer  in  (b)).  Figure 
copyright  2013,  AIP  Publishing,  LLC.  Reproduced  with  permission  from  [18]. 

absorption  related  to  its  associated  ground  state  transition  width  T.  Since  the 
resonances  are  generally  very  narrow  compared  to  the  energy  widths  in  typical 
photon  irradiation  source  spectra,  the  cross  sections  are  typically  given  as  integrated 
cross  sections  (ICS)  in  the  literature  [17,  19]. 

Photoexcitation  by  Bremsstrahlung. 

One  method  of  resonantly  exciting  a  nucleus  into  an  isomer  is  through 
irradiation  by  bremsstrahlung.  Bremsstrahlung  photons  are  continuously 
distributed  in  energy  up  to  the  endpoint  energy  of  the  accelerated  electrons  used  to 
create  them.  If  E  is  photon  energy,  ^  represents  the  spectral  flux  density  of  the 
bremsstrahlung  source,  and  the  number  Nf  of  nuclei  in  a  particular  sample  excited 
into  the  isomer  from  Ni  initial  nuclei  is 


(7) 


Here  the  sum  is  over  all  intermediate  states  that  feed  the  isomer,  identified  by  index 
j.  The  term  ICS fj  is  the  integrated  cross  section  for  population  of  the  isomer  via 
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the  jth  IS  [20], 

This  energy-  and  angle-integrated  cross  section  can  be  evaluated  directly  using 
data  available  in  nuclear  databases  such  as  the  ENSDF.  Dracoulis,  et  al.  identifies 
the  integrated  cross  section  for  an  isolated  (7,7')  resonance  with  photon  wavelength 
A  from  an  initial  state  with  angular  momentum  -Jt  to  an  intermediate  state  with 
angular  momentum  Jj  as 

ICS,,  =  f  a(E)d.E  =  (8) 

E 

with  a  statistical  factor  g  given  by 

2  Jj  +  1 
9  ~  2  Jj  +  1  ‘ 

The  decay  width  of  the  intermediate  state  in  Equation  8  is  E,-,  while  E,-^  is  the 
decay  width  of  all  paths  that  reach  the  isomer  state  and  E is  the  width  for  decay 
from  the  IS  directly  to  the  final  state  [21].  In  the  case  of  photoexcitation,  the  initial 
state  is  the  ground  state  of  the  nucleus  while  the  final  state  is  that  of  the  isomer. 
The  indices  are  reversed  in  the  case  of  isomer  depletion.  The  cross  sections  for 
photoexcitation  and  depletion  are  different  as  a  result  of  the  differing  statistical 
factors  and  the  fact  that  a  given  state  cannot  typically  be  an  IS  for  both  processes. 
Unfortunately  this  formula  has  limited  utility  in  practice,  as  the  natural  decay 
width  T  depends  on  the  level  lifetime  r  according  to  T  =  h/r,  and  these  level 
lifetimes  are  unknown  for  many  levels  in  the  ENSDF.  Instead,  integrated  cross 
sections  can  be  deduced  by  direct  measurement  of  inelastic  photon  excitation,  the 
results  of  which  are  sometimes  published  [17,  20]. 
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2.6  Neutron-Induced  Excitations 


In  the  experiment  described  in  this  work,  (n,  xrij)  reactions  with  x  >  1  were 
used  to  study  the  excited  states  of  rhenium  nuclei.  In  these  reactions,  an  incident 
high-energy  neutron  imparts  energy  to  the  target  nucleus  through  the  formation  of 
a  compound  nucleus.  The  excitation  energy,  which  includes  both  the  kinetic  energy 
and  the  binding  energy  of  the  neutron,  is  distributed  among  the  nucleons  according 
to  equipartition  of  energy.  When  a  nucleon  in  the  compound  nucleus  achieves 
enough  energy  to  tunnel  through  the  potential  barrier  of  the  nucleus,  it  is  emitted 
from  the  nucleus.  In  the  case  where  the  emitted  particle  is  a  neutron,  the  height  of 
the  potential  barrier  is  equal  to  the  neutron  separation  energy  and  this  process  is 
known  as  neutron  evaporation  [22],  For  (rt,  xrr{)  reactions  with  x  >  1,  the 
excitation  energy  of  the  compound  nucleus  is  sufficient  to  overcome  the  separation 
energy  of  x  neutrons  in  the  nucleus,  and  these  neutrons  escape  from  the  nucleus  as  a 
result.  After  neutron  emission,  the  residual  nucleus  is  generally  left  in  an  excited 
state,  and  it  decays  via  y-ray  cascade  or  internal  conversion  to  the  ground  state  or 
to  an  isomer. 

Internal  conversion  is  a  process  that  competes  with  7  decay,  in  which  an  electron 
from  an  inner  atomic  shell  is  ejected  from  the  atom  as  a  means  of  reducing  the 
excitation  energy  of  the  nucleus.  Following  the  emission  of  the  conversion  electron 
from  the  inner  atomic  orbital,  outer  electrons  will  fill  the  vacancy,  emitting  a 
cascade  of  X-rays  or  Auger  electrons  in  the  process  [14].  These  characteristic  X-rays 
appear  prominently  in  most  7-ray  spectra.  The  internal  conversion  coefficient  a  is 
defined  such  that  the  probability  a  given  excited  state  will  decay  via  internal 
conversion  instead  of  7  decay  is  a/(l  +  a).  Conversion  coefficients  are  included  in 
the  calculation  of  partial  7-ray  cross  sections  from  measured  7-ray  yields  for  the 
nuclear  transitions  of  interest  in  this  experiment.  These  conversion  coefficients  are 
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generally  published  in  the  Table  of  Isotopes  [2] ,  or  they  can  be  computed  using  the 
Brlcc  conversion  coefficient  calculator  [23]. 

The  two  primary  reactions  reflected  in  the  data  taken  with  the  GEANIE 
detector  array  were  187Re(n,  n,7)18'Re  and  187Re(n,  2n7)186Re.  Although  (n,xn'y) 
reactions  with  x  >  3  also  occurred  in  the  target,  because  of  the  small  cross  sections 
for  these  reactions  compared  to  those  of  the  x  —  1  and  x  =  2  reactions  there  were 
less  data  available  to  analyze.  Furthermore,  the  focus  of  this  research  in  particular 
was  on  the  structure  of  186Re  revealed  through  the  187Re(n,  2n7)186Re  reaction,  so 
analysis  of  the  data  from  other  reactions  was  a  lesser  priority. 

Reaction  Energetics  and  Cross  Sections. 

Neutron- induced  reactions  have  thresholds  governed  by  the  binding  energy  of 
the  compound  nucleus  and  the  neutron  separation  energy  of  the  target  nucleus.  As 
a  rough  approximation,  the  neutron  separation  energy  of  the  target  nucleus 
represents  a  threshold  below  which  (n,  2 ny)  reactions  are  energetically  forbidden.  In 
this  approximation,  the  187Re  nucleus  must  be  excited  to  an  energy  at  least  as  great 
as  the  neutron  separation  energy  Sn  of  the  18 'Re  nucleus  to  cause  the  reaction  to 
occur  and  create  186Re.  The  neutron  separation  energy  is  given  by  K.  Krane  for  an 
arbitrary  nucleus  z^n  as 

Sn  =  [m  (^Xjv-i)  -  m  (iXN)  +  mn]  c2.  (10) 

For  187Re  this  separation  energy  is  equal  to  8.3  MeV  [14]. 

The  probabilities  of  decay  of  the  compound  nucleus  via  different  reaction 
channels  lead  to  cross  sections  for  the  various  reactions.  The  cross  sections  for  the 
187Re(n,  xn)188_a:Re  reactions  (1  <  x  <  3)  are  plotted  in  Figure  12  using  data 
obtained  from  the  Evaluated  Nuclear  Data  File  (ENDF).  Due  to  the  lack  of 
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experimental  data  for  neutron-induced  reactions  in  187Re,  the  cross  section  data 
adopted  for  the  most  recent  ENDF  database,  ENDF-VII.l,  consists  of  values 
calculated  by  the  TALYS  nuclear  reaction  modeling  code  (further  described  in 
Section  5.3  of  this  document).  These  calculated  cross  sections  are  described  as  a 
good  fit  to  the  available  experimental  data  [24], 


Figure  12.  Cross  sections  for  ( n,n '),  (n,  2n)  and  (n,3n)  reactions  in  18'Re,  plotted 
against  neutron  energy  En.  For  En  >1.0  MeV  the  cross  section  for  radiative  capture 
(n,  7)  is  negligible  compared  to  the  reactions  shown.  Data  for  this  plot  were  obtained 
from  the  ENDF  [24,  25]. 


As  anticipated,  above  the  187Re  neutron  separation  energy  of  8.3  MeV,  the  cross 
section  for  the  18'Re(n,  u')187Re  reaction  decreases  as  the  (n,2n)  reaction  becomes 
possible.  Similarly,  the  cross  section  for  the  187Re(n,  2u)186Re  reaction  decreases  and 
that  of  the  18'Re(n,  3n)185Re  reaction  increases  above  approximately 
2  x  8.3  MeV  =  16.6  MeV.  This  allows  selection  thresholds,  or  cuts,  to  be  defined  at 
neutron  energies  of  ~  10  MeV  and  ~  20  MeV  to  separate  the  reaction  products  by 
reaction  channel  during  the  experimental  data  analysis. 
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2.7  Literature  Review 


The  ENSDF  database  identifies  a  total  of  24  references  for  nuclear  data 
concerning  186Re.  Of  these  24  references,  two  are  compiled  nuclear  data  sheets,  and 
two  are  articles  written  in  German.  Of  the  remaining  original  articles,  two  describe 
the  level  scheme,  and  the  remainder  are  concerned  with  other  nuclear  data  such  as 
the  half  life  or  quadrupole  moments  [7,  15].  In  addition,  since  the  most  recent 
evaluation  in  2002,  two  other  references  have  been  published  that  describe  the  level 
scheme  of  186Re  [26,  27]. 

Nuclear  Structure  and  Isomer  Discovery. 

The  186Re  nucleus  consists  of  75  protons  and  111  neutrons.  By  nature  of  the  fact 
that  it  is  an  odd-odd  nucleus,  its  level  scheme  is  complex,  and  it  has  a  particularly 
high  level  density  at  low  energies.  Regularly  spaced  rotational  and  vibrational 
bands  predicted  by  collective  models  that  are  evident  in  even-even  nuclei  are  not 
immediately  apparent  in  the  186Re  level  scheme.  The  deformed  potential  (Nilsson) 
model  discussed  in  Section  2.2  is  required  to  determine  the  intrinsic  energy  levels  in 
186 Re,  a  prolate  nucleus. 

An  initial  analysis  of  the  nuclear  structure  of  186Re  was  performed  by  Lanier  and 
others  in  1969  [15].  This  was  followed  by  the  discovery  of  a  long-lived  isomer  at 
149(7)  keV  in  1972  by  Secgmiller,  Lindner,  and  Meyer  [7].  This  isomer,  which  has 
spin-parity  Jn  —  8+,  was  observed  to  decay  via  7-ray  emission  to  the  ground  state 
with  a  half  life  of  2  x  105  years.  An  extract  from  the  level  scheme  for  186 Re  is  shown 
in  Figure  13. 

Evident  in  the  figure  is  the  ~  50  keV  transition  by  which  the  ~  150  keV  isomer 
decays  to  the  Jn  =  3“  state  at  99.3  keV.  From  the  99.3  keV  state,  the  nucleus  can 
de-excite  by  direct  7-decay  to  the  ground  state  or  by  cascade  involving  40.3  keV  and 
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Figure  13.  Selected  intrinsic  levels  in  186Re  showing  the  single-nucleon  quantum  num¬ 
bers  in  the  Nilsson  convention  iST7r  [Nnz A]  below  the  horizontal  axis.  The  spherical 
state  spin-parity  J™  from  which  each  state  arises  is  identified  next  to  the  the  respec¬ 
tive  level  in  the  diagram.  Dashed  lines  and  parentheses  identify  tentative  assignments. 
Figure  copyright  1972  Elsevier  B.V.  Reproduced  with  permission  from  [7]. 


59  keV  7-rays.  Transitions  from  the  J71"  =  8+  isomer  to  lower  spin  states  are 
inhibited  by  the  spin  differences  between  these  states. 

A  recent  work  by  Wheldon  and  others  utilized  a  21  MeV  proton  beam  to  excite 
186  Re  nuclei  via  the  transfer  reaction  187Re(p,  d)186Re,  and  identified  30  levels  above 
400  keV  that  had  not  previously  been  observed  [26].  The  new  levels  discovered 
during  the  experiment  were  compared  against  theoretical  models  to  propose  Nilsson 
configurations  for  several  levels,  though  others  are  identified  with  no  proposed 
spin-parity  assignment.  The  data  from  this  experiment  are  included  in  the  XUNDL. 


Recent  Analysis  of  an  Isomer  Feeding  Mechanism. 

A  2014  AFIT  thesis  by  A.  Lerch  examined  spectra  from  an  experiment  involving 
185Re(n,  7)186Re  reactions  performed  at  the  Budapest  Research  Reactor  (BRR).  The 
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results  included  discovery  of  12  levels  and  54  transitions  not  previously  identified, 
bnt  a  transition  feeding  the  149(7)  keV  isomer  level  was  not  among  these  [27].  The 
work  proposed  that  a  possible  mechanism  is  one  in  which  the  isomer  is  fed  by  an  E 2 
transition  from  the  186  keV  level  to  the  149  keV  isomer  level,  but  no  experimental 
or  modeled  data  were  available  to  support  this  hypothesis. 
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III.  Experiment 


3.1  Experimental  Facilities 

The  experiment  took  place  at  the  Los  Alamos  Neutron  Science  Center 
(LANSCE)  Weapons  Neutron  Research  (WNR)  facility  from  10  November  to 
9  December  2014,  utilizing  the  GErmanium  Array  for  Neutron  Induced  Excitations 
(GEANIE)  detector  array. 

LANSCE/WNR. 

The  LANSCE  accelerator  is  a  proton  linear  accelerator  (LINAC)  capable  of 
producing  protons  in  two  distinct  beams  at  energies  up  to  800  MeV.  The  accelerator 
produces  pulses  of  protons  at  a  macropulse  repetition  rate  of  60  or  120  Hz  and 
average  currents  of  0.5  to  4  /xA.  The  beam  is  split  at  locations  along  the  beam  line 
and  pulses  are  delivered  to  various  experiments  at  LANSCE.  Of  the  60  or  120  Hz 
produced  by  the  accelerator  facility,  40  or  100  Hz,  respectively,  is  generally  the 
macropulse  rate  as  it  is  delivered  to  the  GEANIE  target.  The  macropulses,  which 
are  each  625  /is  in  duration,  comprise  numerous  sub-nanosecond  micropulses  spaced 
1.8  /xs  apart. 

For  experiments  at  the  WNR,  the  800  MeV  protons  are  directed  at  one  of  two 
tungsten  targets  to  produce  sources  of  spallation  neutrons  for  the  various  flight 
paths,  illustrated  in  Figure  14.  Each  target  is  located  under  a  large  mound  of 
concrete  that  provides  necessary  shielding. 

Flight  paths  originating  from  target  4  receive  an  unmoderated  flux  of  neutrons 
with  energies  ranging  from  100  keV  to  nearly  600  MeV,  depending  on  the  angle 
between  the  flight  path  and  the  proton  beam  (Figure  15). 
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Figure  14.  The  LANSCE/WNR  facility  consists  of  two  spallation  targets  and  associated 
flight  path  experiments.  GEANIE  uses  neutrons  produced  at  Target  4,  while  the  Blue 
Room  experiments  at  Target  2  involve  exposing  targets  directly  to  the  800  MeV  proton 
beam.  Diagram  authored  by  an  employee  or  employees  of  Los  Alamos  National  Security, 
LLC  (LANS),  operator  of  the  Los  Alamos  National  Laboratory  under  Contract  No.  DE- 
AC52-06NA25396  with  the  U.S.  Department  of  Energy.  Reproduced  with  permission 
from  [28]. 


GEANIE. 

The  GEANIE  spectrometer  is  an  array  of  high-purity  germanium  (HPGe) 
detectors  located  at  60°  right  off  the  proton  beam  line.  The  array  is  used  for  y-ray 
spectroscopy  experiments  that  include  nuclear  structure  investigations  and  cross 
section  measurements  for  various  nuclear  technology  and  national  security 
applications  [28].  A  schematic  of  the  neutron  beam  tube,  fission  chamber  and 
GEANIE  detector  array  is  shown  in  Figure  16. 

Neutron  energies  of  0.6  to  ~  200  MeV  are  available  at  the  target  in  GEANIE 
experiments  [30,  31].  The  neutron  flux  at  the  target  is  measured  using  a  fission 
chamber  positioned  at  the  terminus  of  the  neutron  beam  tube,  18.48  m  from  the 
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Figure  15.  Neutron  flux  spectra  (log-log  scale)  for  target  number  4  at  LANSCE/WNR, 
a  tungsten  spallation  target.  Neutrons  are  produced  with  energies  up  to  800  MeV.  The 
units  of  the  vertical  axis  are  arbitrary  flux  units  defined  such  that  each  flux  distribution 
is  normalized.  Diagram  authored  by  an  employee  or  employees  of  the  Los  Alamos 
National  Security,  LLC  (LANS),  operator  of  Los  Alamos  National  Laboratory  under 
Contract  No.  DE-AC52-06NA25396  with  the  U.S.  Department  of  Energy.  Reproduced 
with  permission  from  [28]. 


spallation  target  [32],  The  fission  chamber  consists  of  an  ionization  chamber  that 
incorporates  iron  foils  impregnated  with  235U  and  238U,  so  that  the  known  fission 
cross  sections  of  these  isotopes  can  be  used  to  calculate  the  neutron  flux  from  the 
counts  of  fission  events  in  the  ionization  chamber  [33].  The  size  of  the  neutron  beam 
arriving  at  the  target  is  adjusted  using  an  iron  collimator,  prior  to  the  beam 
transiting  the  fission  chamber.  For  this  experiment,  the  beam  was  trimmed  to  a 
diameter  of  1/2"  (1.25  cm). 

The  detector  array,  shown  in  Figure  17,  is  positioned  20.34  m  from  the  spallation 
target  and  includes  a  mix  of  planar  and  coaxial  geometry  detectors  to  optimize 
resolution  and  efficiency  across  a  wide  range  of  y-ray  energies.  Coaxial  HPGe 
detectors  have  a  greater  efficiency  at  y-ray  energies  of  ~  400  keV  and  higher,  while 
planar  detectors  have  a  better  low-energy  efficiency.  The  gain  on  the  two  detector 
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Figure  16.  Schematic  diagram  of  the  60°  right  flight  path  from  target  4  at 
LANSCE/WNR  and  the  GEANIE  detector  array.  Diagram  authored  by  an  employee  or 
employees  of  Los  Alamos  National  Security,  LLC  (LANS),  operator  of  the  Los  Alamos 
National  Laboratory  under  Contract  No.  DE-AC52-06NA25396  with  the  U.S.  Depart¬ 
ment  of  Energy.  Reproduced  with  permission  from  [29]. 


types  is  set  in  the  hardware  to  reflect  this  difference.  The  planar  gain  is  set  to 
achieve  full-scale  (8192  channel)  deflection  at  E1  —  1.0  MeV,  while  the  coaxial  gain 
is  set  such  that  full-scale  deflection  occurs  at  E1  =  4.0  MeV.  Because  of  neutron 
damage  to  the  semiconductor  crystals,  some  detectors  failed  to  produce  usable 
spectra,  which  impacted  the  quantity  of  7  —  7  coincidence  data  produced. 
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Figure  17.  The  GEANIE  detector  array  at  LANSCE/WNR,  with  the  enriched  ls'Re 
sample  mounted  in  the  center  of  the  array.  The  beam  originates  from  the  far  right  of 
the  array  in  the  photograph.  This  picture  was  taken  with  the  detector  array  opened, 
so  approximately  2/3  of  the  18  installed  detectors  are  visible. 


The  detector  array  has  mounts  for  26  HPGe  detectors,  which  provide  nearly  47t 
coverage  of  the  target.  For  the  rhenium  experiment,  these  detectors  were  mounted 
in  the  array  according  to  the  diagram  shown  in  Figure  18.  Two  positions  in  the  20 
circumferential  locations  were  left  vacant  in  this  experiment  due  to  detector 
non-availability.  Of  the  18  detectors  installed  during  this  experiment,  all  were 
Compton-suppressed  using  bismuth  germanium  oxide  (BGO)  detectors.  The  6 
overhead  detector  mounts,  which  can  accommodate  unsuppressed  HPGe  detectors, 
were  also  left  empty  for  the  rhenium  experiment  due  to  detector  non-availability. 
For  this  particular  experiment,  the  detector  array  consisted  of  18  serviceable 
detectors,  a  mix  of  8  planar  and  10  coaxial  detectors. 
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BEAM 


(a)  Top-down  view  (b)  Side  view 

Figure  18.  Diagram  of  the  HPGe  detector  locations  in  the  GEANIE  array  (not  to  scale), 
viewing  from  (a)  top  down  perspective  and  (b)  side  perspective,  with  the  neutron  beam 
coming  out  of  the  page.  The  diagrams  identify  the  18  detectors  that  were  installed  and 
functional  in  the  array  at  the  time  of  the  experiment.  Letters  identifying  the  detector 
are  referenced  in  Table  1  with  the  azimuthal  angle  6  and  angle  of  elevation  <p. 


The  polar  and  azimuthal  detector  locations  relative  to  the  target  are  listed  in 
Table  1  [29].  Remarks  in  the  table  identify  problems  with  certain  detectors 
encountered  during  the  experiment  and  data  analysis,  described  further  in 
Chapter  IV. 
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Table  1.  GEANIE  detector  positions  by  azimuthal  angle  0  and  angle  of  elevation  <fi 
relative  to  beam  axis.  The  mounts  for  detectors  12,  13,  and  21  26  were  empty  during 

this  experiment. 


Detector  Type 

0 

9 

Distance 

Remarks 

1 

Planar 

29.0° 

-152.8° 

14.415 

cm 

2 

Planar 

-29.0° 

-154.0° 

14.442 

cm 

3 

Planar 

29.0° 

1.0° 

14.379 

cm 

4 

Coaxial 

-29.0° 

1.2° 

14.773 

cm 

5 

Planar 

0.0° 

26.5° 

14.308 

cm 

Double  peaks  in  TDC  spectrum 

6 

Coaxial 

29.0° 

53.0° 

14.455 

cm 

7 

Coaxial 

-29.0° 

53.5° 

14.435 

cm 

8 

Planar 

0.0° 

78.5° 

14.917 

cm 

9 

Planar 

29.0° 

102.0° 

14.237 

cm 

10 

Planar 

-29.0° 

102.5° 

14.288 

cm 

11 

Coaxial 

0.0° 

129.5° 

14.237 

cm 

No  counts  after  run  9773 

14 

Coaxial 

0.0° 

-25.2° 

14.392 

cm 

Double  peaks  in  TDC  spectrum 

15 

Coaxial 

29.0° 

-51.1° 

14.392 

cm 

Double  peaks  in  TDC  spectrum 

16 

Coaxial 

-29.0° 

-51.0° 

13.846 

cm 

Unstable  ADC  gain 

17 

Coaxial 

0.0° 

-76.9° 

14.442 

cm 

18 

Coaxial 

29.0° 

-101.7° 

14.176 

cm 

19 

Coaxial 

-29.0° 

-102.0° 

14.308 

cm 

20 

Planar 

0.0° 

-128.0° 

14.161 

cm 

Compton  Suppression. 

The  response  function  of  HPGe  detectors,  like  all  7-ray  detectors  of  finite  size, 
includes  a  continuum  of  partial  energies  deposited  in  the  detector  corresponding  to 
the  energies  of  Compton  electrons  scattered  by  the  incident  7-ray.  The  Compton 
continuum  extends  up  to  energies  of 


EP-  = 


2  E% 


2  E1  +  mec 2  ’ 


(a) 
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where  EL,  is  the  energy  of  the  7-ray,  me  is  the  mass  of  the  electron,  and  c  is  the 
speed  of  light.  Equation  11  identifies  the  maximum  energy  attainable  by  a  Compton 
electron,  when  the  scattering  angle  is  9  =  it  [14].  In  7-ray  spectroscopy,  the  presence 
of  the  Compton  continna  from  higher  energy  7-rays  represents  a  source  of 
background  that  might  obscure  the  full  energy  photopeaks  from  lower  energy  7-rays. 

A  Compton  suppression  system  was  designed  to  alleviate  this  problem.  In  the 
design  used  in  the  GEANIE  detectors,  a  bismuth  germanium  oxide  (BGO)  detector 
provides  a  shield  surrounding  each  HPGe  detector.  BGO  detectors  are  selected  for 
this  role  because  of  their  high  relative  detection  efficiency,  due  to  the  presence  of 
liigh-Z  bismuth  in  the  detector  crystals.  If  a  Compton-scattered  photon  is  measured 
in  the  BGO  detector  at  the  same  time  that  a  Compton  electron  is  measured  in  the 
HPGe  detector,  then  the  signal  is  rejected.  This  anticoincidence  setting  ensures  that 
the  only  signals  sent  from  the  suppressed  HPGe  detectors  to  the  ADCs  are  full 
photopeaks. 

3.2  Target  Sample  Preparation 

Two  targets  of  approximately  1  gram  each  were  used  in  this  experiment,  a 
rhenium  sample  enriched  in  18'Re  and  a  sample  of  natural  rhenium.  Both  were 
ordered  by  Ecopulse,  Inc.  under  a  contract  maintained  by  the  Army  Research 
Laboratory  (ARL).  The  rhenium  metal  powder  samples  were  hot  vacuum  pressed  at 
Oak  Ridge  National  Laboratory  (ORNL)  and  sent  directly  to  LANSCE  in  October 
and  November,  2014. 

Enriched  18  Re  Sample. 

The  enriched  rhenium  sample  was  purchased  from  ORNL,  where  it  was 
isotopically  enriched  to  99.52%  18'Re.  The  encapsulated  sample  consisted  of 
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981.7  mg  rhenium  metal  powder. 


nat Re  Sample. 

A  natural  rhenium  (37.4%  185 Re  and  62.6%  18'Re)  sample  was  used  in  the  last 
day  of  the  experiment  and  consisted  of  1036.2  mg  of  metal  powder.  This  sample  was 
purchased  from  Alfa  Aesar  and  shipped  to  ORNL  for  encapsulation. 

Target  Holders. 

Polycarbonate  target  holders  were  3-D  printed  in  two  pieces  at  ARL,  as  shown 
in  Figure  19.  These  two  pieces  were  held  together  by  setscrew,  and  contained  the 
vacuum-pressed  metal  pill  samples  after  they  were  fashioned. 


Source  Holder  BASE 


Source  Holder  CAP 


Figure  19.  Schematic  of  the  3D-printed  holder  used  to  contain  the  pills  of  pressed 
rhenium  metal  powder  during  irradiation. 


The  target  holder  was  made  to  eliminate  the  possibility  of  leakage  of  the 
rhenium  powder  during  the  experiment,  as  trace  amounts  of  radioactive  186 Re  would 
be  created  by  (n,  2 n)  reactions.  However,  concerns  of  leakage  from  the 
encapsulation  led  to  special  safety  procedures  mandated  by  Radiation  Control 
Technicians  (RCTs)  at  LANL,  in  which  the  sample  was  required  to  be  tested  for  the 
presence  of  the  /3-emitting  isotope  186Re  each  time  the  target  sample  was  removed 
from  the  detector  array. 
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3.3  Experiment  Summary 


The  21  days  of  beam  time  originally  requested  for  the  experiment  was  approved 
in  two  parts:  the  first  part  from  10  -  23  November  2014  and  the  remainder  from 
4-8  December  2014.  Due  to  the  timing  of  other  experiments,  beam  was  available 
from  3-9  December  during  the  second  window.  During  the  approved  beam  time 
windows,  two  related  experiments  were  scheduled  involving  the  GEANIE  detector 
array.  The  first  involved  the  rhenium  sample  enriched  in  187Re  and  the  other 
involved  the  natural  rhenium  sample. 

Enriched  18  Re  Target. 

Using  the  enriched  18'Re  target,  7-ray  spectra  from  the  187Re(n,  n/7)187Re  and 
187Re(n,  2n7)186Re  reactions  were  obtained  over  20  days,  encompassing  the 
10-23  November  window  and  the  first  6  days  of  the  3-9  December  window.  The 
experiment  was  initially  planned  anticipating  a  100  Hz  macropulse  repetition  rate 
produced  by  the  LANSCE  accelerator,  but  beam  maintenance  issues  led  to  a  lower 
repetition  rate  of  40  Hz  and  an  average  proton  current  of  approximately  1200  11A 
from  10-23  November.  Additionally,  a  bad  filament  in  the  accelerator  led  to  an 
unplanned  28-hour  shutdown  from  17-18  November.  As  a  result  of  the  lower  current 
and  the  unplanned  outage,  the  decision  was  made  to  extend  the  experiment 
involving  the  enriched  target  to  16  days  from  the  8  days  originally  planned. 

During  the  planned  24  November-2  December  shutdown,  accelerator 
maintenance  personnel  increased  the  macropulse  repetition  rate  to  100  Hz,  and  the 
experiment  resumed  on  3  December  using  the  enriched  target.  From 
3-9  December,  the  accelerator  operated  continuously  at  100  Hz  and  an  average 
proton  beam  current  of  ~  3500  11A. 
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Normalization  Runs  with  56Fe. 


From  6-8  December,  the  enriched  rhenium  sample  was  sandwiched  between 
two  0.127-nnn  thick  na*Fe  foils  and  irradiation  continued.  The  data  obtained  from 
these  runs  allows  measurement  of  the  7-ray  yield  for  the  846.8  keV  transition  in  56Fe 
from  the  56Fe(n,  n/7)56Fe  reaction,  relative  to  the  7-ray  yields  from  neutron-induced 
reactions  in  rhenium.  The  well-known  partial  7-ray  cross  section  for  the  846.8  keV 
transition  in  56Fe  is  used  to  normalize  the  partial  7-ray  cross  sections  generated  for 
the  rhenium  isotopes,  so  that  the  rhenium  cross  sections  can  be  compared  with 
those  produced  by  theoretical  models. 

nat Re  Target. 

Using  the  sample  of  naiRe,  7-ray  spectra  from  185Re(n,  n/7)185Re  and 
185 Re(n,  2/17)  184Re  reactions  were  obtained  in  addition  to  the  reactions  involving 
18 'Re  over  a  period  of  24  hours  from  8-9  December.  The  intent  of  irradiating  the 
natural  rhenium  target  was  to  generate  7-ray  excitation  functions  from  the 
185 Re(n,  2/17)  184Re  channel  for  comparison  with  the  excitation  functions  obtained 
from  similar  transitions  in  18'Re(n,  2n7)186Re,  which  will  be  described  in  more 
detail  in  Section  5.4. 

Source  Data. 

When  beam  operations  were  stopped,  calibration  sources  were  placed  in  the 
detector  array  for  the  purpose  of  obtaining  spectra  for  energy  and  efficiency 
calibration  of  the  detectors.  A  total  of  9  calibration  sources  were  used  to  cover  a 
wide  range  of  7-ray  energies:  152Eu,  57Co,  133Ba,  60Co,  22Na,  109Cd,  137Cs,  54Mn  and 
241  Am.  Data  were  collected  in  20-minute  runs  for  each  of  these  sources  during  the 
unscheduled  stoppage  on  17  November.  Additional  source  data  using  the  152Eu  and 
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241  Am  sources  were  also  collected  during  the  scheduled  24  November  -  3  December 
shutdown. 

The  data  from  these  sources,  with  the  152Eu  source  as  the  primary  reference  for 
energy  calibration,  were  analyzed  separately  to  produce  energy  and  efficiency 
calibration  hies  for  use  in  the  RadWare  programs  gf3  and  escl8r,  described  in 
further  detail  in  Chapter  IV. 

Detector  Performance. 

During  the  experiment,  14  of  the  18  installed  HPGe  detectors  (7  planar  and 
7  coaxial)  produced  usable  spectra  that  were  included  in  the  data  analysis.  The 
remaining  4  detectors  developed  problems  that  prevented  the  inclusion  of  the 
spectra  they  produced  in  the  final  data  set.  Those  issues  are  identified  in  the 
remarks  column  of  Table  1,  and  the  reasons  for  their  non-inclusion  are  as  follows: 

•  One  of  the  installed  planar  detectors  (5)  and  2  coaxial  detectors  (14,  15)  were 
identified  as  producing  time  spectra  with  anomalous  double  peaks.  This  effect 
prohibits  the  correlation  of  neutron  energy  En  to  the  TDC  channel,  so  these 

3  detectors  were  eliminated  from  the  set  during  analysis  of  the  data. 

•  One  coaxial  detector  (16)  had  a  highly  unstable  gain  that  caused  significant 
problems  during  gain  matching,  so  it  was  eliminated  from  the  final  detector 
set. 

Of  the  remaining  detectors,  one  coaxial  detector  (11)  ceased  producing  counts 
altogether  1  week  after  the  experiment  began.  Based  on  this  detector’s  location 
opposite  the  neutron  source  it  is  likely  that  neutron  damage  resulted  in  the 
detector’s  early  failure.  Prior  to  the  failure,  however,  the  detector  produced  usable 
spectra  which  were  included  in  the  final  data  set. 
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There  were  11  single-channel  wide  lines  in  ADCs  1-6,  which  ultimately  included 
3  planar  detectors  and  2  coaxial  detectors,  that  were  attributed  to  the  hardware,  as 
all  these  detectors  were  connected  to  the  same  data  bus.  These  spurious  lines  were 
noted  for  reference  during  the  data  analysis  to  remove  the  possibility  of 
mis-identifying  them  as  7-ray  peaks. 

During  the  3-8  December  runs,  Compton  suppression  was  lost  to  all  the 
detectors  in  the  array  for  a  period  of  3  days.  The  problem  was  eventually  traced  to 
a  tripped  power  supply,  which  affected  the  entire  bus  from  which  the  BGO 
suppression  shields  received  their  power.  The  experiment  was  completed  by  the 
time  the  problem  was  discovered,  so  the  7-ray  spectra  obtained  from  the  natRe 
sample  included  a  Compton  background  not  present  in  the  earlier  runs. 
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IV.  Results  and  Analysis 


The  experiment  was  unique  in  two  ways:  it  was  the  first  reported  experiment  to 
study  18 'Re  with  an  (n,n')  reaction,  and  the  first  to  study  186 Re  with  a  (n,2n) 
reaction.  As  a  result,  it  was  expected  that  the  experiment  would  result  in  the 
discovery  of  new  7-ray  transitions  in  both  isotopes.  The  initial  goals  of  data 
analysis  following  the  experiment  were  to  identify  7-ray  peaks  in  the  ADC  spectra 
produced  by  the  HPGe  detectors  attributable  to  new  transitions  in  186 Re  and  18'Re. 


4.1  Data  Acquisition  and  Online  Analysis 


Data  acquired  by  GEANIE  were  first  analyzed  online,  while  the  experiment  was 
running,  prior  to  being  saved  for  offline  data  analysis.  The  steps  involved  in  this 
process  are  shown  schematically  in  Figure  20. 


o 
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Figure  20.  Flow  chart  showing  the  process  of  data  acquisition,  storage  and  transfer. 
The  necessary  files  for  offline  data  analysis  consisted  of  ROOT  histograms  (.root)  and 
MIDAS  raw  data  (.mid.gz)  and  online  database  (.odb)  information. 
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Data  Acquisition  Framework. 


In  the  online  phase  of  the  analysis,  7-ray  energy  and  neutron  time-of-flight 
(TOF)  data  were  acquired  by  analog-to-digital  converters  (ADCs)  and 
time-to-digital  converters  (TDCs).  The  ADC  and  TDC  signals  were  sent  to  Versa 
Module  Europa  (VME)  and  Computer  Automated  Measurement  and  Control 
(CAMAC)  servers  that  interfaced  with  a  Linux  workstation  operating  as  a  Data 
Acquisition  (DAQ)  server.  The  VME  server  handled  data  sent  to  the  VME  bus, 
which  included  ADC  and  TDC  signals  from  each  of  the  HPGe  detectors  in  the 
GEANIE  array  and  the  fission  chambers.  Signals  from  scalers,  which  recorded  clock 
timing  data  and  the  total  ADC  counts  over  each  10-second  period  of  runtime 
(necessary  for  calculating  detector  deadtime),  were  sent  through  the  CAMAC  server. 

All  experiments  at  LANSCE/WNR,  including  GEANIE,  use  the  Maximum 
Integrated  Data  Acquisition  System  (MIDAS)  software  package  as  the  DAQ 
framework  to  interface  with  the  VME  and  CAMAC  hardware.  MIDAS  was 
developed  at  the  Paul  Scherrer  Institute,  Switzerland  for  particle  physics 
experiments,  and  is  open-source  software. 

MIDAS  Analyzer. 

Data  from  each  7-ray  event  were  compiled  on  the  MIDAS  server  using  an 
analyzer  program  written  at  LANSCE,  which  performed  online  analysis  of  the  VME 
and  CAMAC  signals  to  produce  event  data  that  makes  up  each  run.  After  each  run 
was  completed,  the  MIDAS  analyzer  packaged  and  compressed  the  run  hies,  and 
data  from  each  run  were  saved  to  disk.  The  hies  produced  at  this  point  included  the 
raw  data  from  each  7-ray  event  in  the  run,  saved  in  a  proprietary  MIDAS  format 
and  compressed  with  the  hie  extension  .mid.gz,  histograms  that  display  the  various 
ADC  and  TDC  spectra  in  a  ROOT  format,  and  an  online  database  ( .  odb)  hie 
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containing  various  header  information  about  the  run  for  use  by  the  MIDAS  analyzer 
program.  The  largest  of  these  hies  from  each  run  was  generally  the  .mid.gz  hie, 
which  ranged  in  size  for  each  run  from  several  hundred  megabytes  to  over  a 
gigabyte,  depending  on  the  run  duration. 

The  .mid.gz,  .odb  and  .root  hies  were  saved  locally  to  disk  on  the  MIDAS 
server  automatically  following  each  run,  and  the  next  run  started  automatically 
according  to  the  preset  run  times.  The  run  times  were  adjusted  throughout  the 
experiment  to  minimize  data  loss  while  also  limiting  the  total  number  of  data  hies 
created.  The  online  MIDAS  analyzer  experienced  frequent  crashes,  particularly 
when  the  data  rate  was  high,  as  was  the  case  during  beam  operations  at  the  higher 
100  Hz  macropulse  repetition  rate.  Setting  the  run  time  to  a  smaller  value  insured 
against  data  loss  when  the  analyzer  crashed,  as  the  program  would  stop  and  restart 
itself  automatically  at  the  prescribed  time  and  clear  the  buffer  in  the  process.  In 
general,  run  times  varied  between  20  minutes  and  6  hours. 

Data  Transfer  and  Storage. 

Periodically  throughout  the  experiment,  the  hies  were  uploaded  to  the  LANL 
public  FTP  server,  from  which  they  were  downloaded  to  an  external  hard  drive.  At 
the  end  of  the  experiment,  the  total  data  accumulated  on  the  hard  drive  was 
62.4  GB  in  size,  which  included  data  from  441  runs. 

4.2  Offline  Data  Analysis  Methodology 

The  MIDAS  event  hies  were  converted  into  y-ray  energy  (Ey)  versus  neutron 
energy  (En)  matrices  during  the  offline  phase  of  data  analysis.  These  E1  vs.  En 
matrices  were  then  parsed  into  y-ray  spectra  gated  on  specihc  neutron  energy 
ranges. 
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Software  Alternatives  for  Data  Analysis. 


In  order  to  convert  MIDAS  event  files  to  viewable  data  structures,  they  had  to  be 
converted  to  be  viewed  by  an  appropriate  offline  analysis  package.  Two  alternatives 
were  available  for  this  analysis:  ROOT  and  TScan.  With  TScan,  the  7-ray  spectra 
are  analyzed  using  programs  included  in  the  RadWare  software  package. 

ROOT. 

ROOT  is  an  open-source  collection  of  C++  libraries  developed  and  maintained 
at  CERN.  The  well-documented  libraries  include  statistics  and  peak  fitting  tools 
that  allow  users  to  analyze  spectra  directly.  It  is  considered  a  standard  for  data 
analysis  in  many  nuclear  and  particle  physics  experimental  facilities.  The  MIDAS 
analyzer  is  capable  of  producing  ROOT  trees  and  histograms  in  an  offline  mode,  so 
ROOT  can  be  used  to  perform  spectral  analysis  of  GEANIE  data  [34,  35].  Of 
interest  in  this  experiment,  however,  was  the  generation  of  7  —  7  coincidence 
matrices,  for  which  there  were  no  tools  built  in  ROOT. 

TScan. 

The  TScan  software  package  was  written  by  H.  Jin  at  Rutgers  University  to 
perform  tape  scanning  and  manipulate  matrices  [36].  It  has  been  used  for  data 
analysis  in  numerous  GEANIE  experiments  [1,  29,  32,  37,  38,  39].  Though  TScan  is 
no  longer  maintained,  it  can  still  be  compiled  and  run  on  a  32-bit  Linux  operating 
system.  Numerous  subroutines  developed  specifically  for  analyzing  data  from 
GEANIE  experiments  have  been  written  by  W.  Younes  of  Lawrence  Livermore 
National  Laboratory,  including  a  tool  for  generating  coincidence  matrices  [40].  As  a 
result,  TScan  was  adopted  for  data  analysis  in  this  experiment  despite  its  legacy 
status. 
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Rad  Ware. 


RadWare  is  an  open-source  software  package  for  analyzing  one-dimensional 
(singles)  7-ray  spectra  and  constructing  nuclear  level  schemes  from  7  —  7 
coincidence  data,  among  other  functions  [41].  It  was  developed  by  D.  Radford, 
currently  of  ORNL,  and  it  is  used  widely  throughout  the  nuclear  physics 
community.  The  spectra  and  coincidence  matrices  produced  by  TScan  and 
associated  subroutines  were  analyzed  using  the  gf3  and  escl8r  programs  included 
with  the  current  RadWare  version  rw05. 

Data  Analysis  Using  TScan  and  RadWare. 

After  the  experiment  was  completed  and  the  data  hies  downloaded  from  the 
FTP  server,  the  run  hies  were  converted  from  the  MIDAS  format  to  a  hie  format 
readable  by  TScan.  After  this  conversion  was  complete,  the  ADC  spectra  from  the 
14  detectors  were  summed  and  gain-matched  using  TScan  to  ensure  peaks  from  all 
detectors  were  properly  aligned  before  singles  spectra  and  7  —  7  coincidence 
matrices  could  be  further  analyzed.  The  steps  involved  in  this  process  are 
illustrated  in  the  how  chart  shown  in  Figure  21. 
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Figure  21.  Flow  chart  describing  the  steps  involved  in  the  offline  analysis  of  data  from 
the  GEANIE  experiment.  Grey  boxes  identify  the  input  files,  which  contain  data  or 
parameters,  used  by  the  programs  (shown  in  orange). 

Run  file  conversion  to  TScan  format. 

After  the  compressed  MIDAS  run  files  in  the  .  mid .  gz  format  were  unzipped, 
they  were  dumped  into  a  hexadecimal  file  using  the  MIDAS  subroutine  ndump.  This 
hexadecimal  file  was  then  sent  as  input  to  the  Perl  scripts  goode  and  itob,  which 
converted  the  data  into  event  files  ( .  evt)  readable  by  TScan.  A  BASH  shell  script 
was  written  to  automate  this  process,  which  generally  required  1-2  minutes  to 
process  each  event  file. 

TScan  presort  routine. 

The  TScan  program  was  compiled  from  a  presort  routine  sort .  c,  which 
contained  the  sorting  parameters  used  to  sum  the  data  from  each  of  the  detectors 
into  a  matrix.  The  result  of  compiling  the  presort  routine  sort .  c  is  the  TScan 
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executable,  which  summed  the  numerous  event  hies  to  produce  a  single 
2-dimensional  E1  vs.  En  matrix  in  the  .  gmt  hie  format.  The  single  .  gmt  hie 
contained  all  the  spectral  data  from  the  experiment,  formatted  per  the  specifications 
defined  in  the  presort  routine. 

The  conversion  from  neutron  TOF  to  energy  was  included  in  the  presort  routine. 
This  step  involved  application  of  the  relativistic  formula 


F  _  E0 

"  y/1  -U2/C2’ 


(12) 


where  Eq  =  939.57  MeV  is  the  neutron  rest  energy,  v  is  the  neutron  velocity  and  c  is 
the  speed  of  light  in  a  vacuum  [14,  39].  Since  the  distance  between  the  proton  beam 
pick-off  to  the  detector  array  was  known  to  be  20.34  m,  the  neutron  energy  En 
could  be  deduced  from  the  TOF.  Also  contained  in  the  presort  routine  were  the 
parameters  necessary  to  align  the  TDC  spectra  from  each  detector  to  the  same 
channel,  so  that  this  TOF-to-F„  conversion  was  consistent  across  all  the  detectors. 

TOF  measurement  involved  a  TDC  measurement  of  the  time  between  the  start 
signal,  generated  by  the  proton  beam  pick-off,  and  the  stop  signal,  which  occurred 
upon  the  7  flash  in  the  HPGe  detectors  from  each  micropulse.  The  TDC  spectrum 
structure  is  shown  in  Figure  22(a).  The  TDC  spectra  from  each  detector  exhibited  a 
447-channel  periodicity  according  to  the  beam’s  constant  micropulse  spacing,  which 
was  used  in  the  presort  routine  to  calculate  the  neutron  time  of  flight  (TOF)  that 
accompanies  each  7-ray  signal.  The  output  from  this  step  is  the  TOF  spectrum 
shown  in  Figure  22(b). 

The  neutron  energies  En  —  1,  10  and  25  MeV  highlighted  in  the  inset  of 
Figure  22(b)  were  calculated  from  the  TOF  using  Equation  12. 
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(a)  TDC  spectrum 


Figure  22.  Time  spectra  generated  by  the  planar  detectors.  The  time  spectra  from  the 
coaxial  detectors  are  similar  to  those  shown  for  the  planar  detectors.  Panel  (a)  shows 
the  raw  TDC  output  from  detector  1.  The  constant  micropulse  spacing  in  the  TDC 
spectrum  was  used  to  collapse  the  TDC  spectra  into  the  neutron  time-of-flight  spectrum 
shown  in  panel  (b),  which  includes  the  TDC  data  from  all  the  planar  detectors.  The 
neutron  energies  at  which  the  gates  were  applied  in  the  data  analysis  are  represented 
by  the  vertical  lines  in  the  inset.  The  single-channel  wide  vertical  spike  at  channel 
1080  is  a  spurious  signal  generated  by  the  TDC  hardware,  and  must  be  corrected  when 
calculating  dead  time  corrections. 


Gain  matching. 

The  ADCs  output  7-ray  energies  (f?7)  in  terms  of  counts  per  channel,  and  had 
to  be  gain-matched  to  ensure  peaks  from  each  detector  were  aligned  at  the  same 
channel.  Due  to  the  different  hardware  gain  settings  for  the  planar  and  coaxial 
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detectors,  the  two  detector  types  were  gain-matched  separately.  The  gain  matching 
procedure  was  complicated  by  the  fact  that  over  the  month-long  experiment  the 
gain  for  certain  detectors  drifted  by  significant  amounts. 

The  gain  matching  was  performed  iteratively  using  small  subsets  of  the  data  to 
insure  against  significant  gain  and  baseline  drifts  that  might  result  in  unacceptable 
peak  broadening.  After  summing  the  events  from  20-30  event  hies  into  a  single  .  gmt 
matrix  hie,  the  matrix  was  parsed  using  the  RGMT  (read  .gmt)  subroutine  from 
TScan  into  one-dimensional  7-ray  spectrum  ( .  spe)  hies.  These  spectra  were 
examined  for  each  of  the  14  detectors  using  the  gf3  program  from  the  RadWare 
package.  The  peak  locations  by  channel  were  recorded  and  used  to  generate 
corrected  gain  and  baseline  parameters  according  to  the  following  formulas  [40]: 

Gain  =  (13) 

x2  -  Xi 

Shift  =  i/2  —  Gain  x  x2 . 

Here  X\  and  x2  are  the  channel  locations  of  two  peaks  in  the  spectrum  with  known 
energy,  while  ij\  and  y2  are  the  desired  channel  locations  of  those  peaks  in  the 
reference  spectrum.  For  the  gain  matching  of  the  planar  detectors,  the  reference 
peaks  selected  were  the  61.14  keV  Kq,!  X-ray  from  rhenium  (see  Table  6)  and  the 
511  keV  annihilation  peak,  which  were  assigned  reference  locations  based  on  a  gain 
of  8.192  channels/keV.  The  peaks  used  for  gain  matching  the  coaxial  detectors  were 
the  134.247  keV  7-ray  from  187Re  and  the  511  keV  e+e_  annihilation  peak,  using  a 
reference  gain  of  2.048  channels/keV.  The  corrected  parameters  for  each  subset  of 
the  data  were  provided  as  input  to  TScan  using  the  hies  gain .  cmd  and  shift .  cmd 
before  summing  the  subsets  into  a  single  master  .  gmt  matrix. 


51 


Singles  spectra  and  coincidence  matrices. 


The  RGMT  program,  as  described  above,  generated  one-dimensional  7-ray  spectra 
from  the  final  .  gmt  matrix.  The  spectra  generated  in  this  way  included  individual 
ADC  spectra  for  each  of  the  14  detectors,  as  well  as  summed  spectra  incorporating 
the  counts  from  all  of  the  planar  or  coaxial  detectors.  The  summed  spectra, 
provided  the  gain  matching  was  performed  properly,  resulted  in  the  best  counting 
statistics  and  thus  the  best  opportunity  for  observation  of  new  7-ray  transitions. 
These  summed  spectra,  gated  on  selected  neutron  energy  ranges  as  described  in 
Section  4.2,  were  the  focus  of  the  data  analysis. 

Included  in  the  TScan  program  is  an  executable  gmt2mat  for  extracting  from  the 
.gmt  matrix  7  —  7  coincidence  matrices  in  the  .mat  format  preferred  by  RadWare. 
The  RadWare  distribution  includes  a  program  escl8r  that  is  designed  to  examine 
these  2-dimensional  coincidence  matrices  and  produce  level  schemes  from  the 
coincidence  data. 

Neutron  Energy  Gates. 

Gates  on  the  neutron  energies  were  specified  in  the  TScan  presort  routine 
sort .  c  to  generate  ADC  spectra  showing  7-rays  that  originated  primarily  from 
(n,  n'y)  or  (n,  2 717)  reactions,  according  to  the  thresholds  of  the  cross  sections 
plotted  in  Figure  12.  The  neutron  energy  gates  used  to  select  7-rays  that  originated 
from  the  187Re(n,  n/7)187Re  reaction  channel  were  1.0  <  En  <  10  MeV.  For  the 
187Re(n,  2n7)186Re  reaction  channel,  Figure  12  identifies  that  this  reaction 
dominates  the  (n,  n'y)  and  (n,  3777)  reactions  in  the  energy  range 
10  <  En  <  18  MeV.  Gates  of  10  <  En  <  25  MeV  were  chosen  to  increase  the 
number  of  counts  in  each  channel,  and  known  7-ray  peaks  due  to  185Re  were 
subtracted  out  to  leave  only  contributions  from  186Re. 


52 


The  1.8  ns  delay  between  micropulses  is  a  sufficiently  long  time  for  the 
high-energy  neutrons  from  one  of  the  sub-nanosecond  micropulses  to  overtake  the 
low-energy  neutrons  from  a  previous  micropulse,  an  effect  known  as  frame  overlap 
[32],  This  must  be  accounted  for  in  the  data  analysis  by  selection  of  an  appropriate 
neutron  low-energy  cutoff  of  650  keV  or  greater.  For  this  experiment,  the 
low-energy  cutoff  chosen  was  1.0  MeV,  which  was  sufficient  to  select  out  the  slow 
neutrons  that  would  be  subject  to  frame  overlap  and  remove  any  ambiguity  in  the 
neutron  energy.  This  low-energy  cutoff  also  had  the  effect  of  eliminating  y-rays 
produced  in  the  187Re(n,  y)188Re  reaction  channel,  because  the  (n,  y)  cross  section  is 
negligible  compared  to  that  of  the  (n,  n')  reaction  at  neutron  energies  above  1  MeV. 

With  the  neutron  energy  gates  set  as  described  above,  the  neutron  energy  bins 
in  the  E1  vs.  En  matrix  were  projected  onto  the  E^  axis  to  generate  summed 
projection  spectra  consisting  only  of  7-rays  correlated  to  neutrons  in  the  desired 
energy  ranges.  These  projection  spectra  are  shown  in  Figures  23  and  24  for  the 
planar  and  coaxial  detectors,  respectively. 

Energy  Calibration. 

Energy  calibration  was  performed  using  the  Source  and  ENCAL  executables  from 
the  RadWare  software  package.  The  primary  calibration  source  for  both  planar  and 
coaxial  detectors  was  152Eu.  Because  of  the  different  gain  settings  for  the  two 
detectors,  different  7-ray  peaks  were  used  for  the  two  calibrations:  7  peaks  with 
energies  between  121  keV  and  867  keV  were  used  for  the  planar  detector  calibration, 
while  10  peaks  with  energies  ranging  from  121  keV  to  1408  keV  were  used  to 
calibrate  the  coaxial  detector  spectrum.  The  ENCAL  subroutine  produced  planar  and 
coaxial  calibration  files  based  on  linear  least-squares  fits  to  the  known-energy  peaks 
in  152Eu.  Further  analysis  is  required  to  quantify  the  errors  inherent  in  the  energy 
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calibration  process,  although  they  were  estimated  to  be  on  the  order  of  0.1  keV. 


Peak  Fitting  in  RadWare  (gf3). 

Peaks  in  the  projection  spectra  were  fitted  using  the  least-squares  peak  fitting 
algorithm  in  gf3.  The  gf3  program  is  specifically  designed  for  analysis  of  7-ray 
spectra  produced  by  HPGe  detectors,  and  it  fits  peaks  according  to  the  unique 
response  function  of  an  HPGe  detector  [41].  For  this  experiment,  each  peak  was 
fitted  with  a  pure  Gaussian  and  peak  centroid  errors  due  to  the  fitting  process  were 
quantified.  These  errors  are  reflected  in  Tables  2  through  5.  Counts  under  each 
peak  were  determined  by  subtracting  the  background,  estimated  over  a  range  of 
channels  bounding  the  peak. 

4.3  Projection  Spectrum  Analysis 

Of  primary  interest  in  this  experiment  was  evidence  of  transitions  feeding  the 
149(7)  keV  isomer  in  186Re,  though  new  transitions  in  1S7Re  were  also  sought.  Peaks 
in  the  projection  spectra,  gated  on  the  neutron  energies  described  in  Section  4.2, 
were  fitted  using  gf3,  and  the  energies  of  these  peaks  were  compared  with  known 
7-ray  transitions  listed  in  the  ENSDF  [10]. 
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Planar  Detector  Projection  Spectra. 


The  projection  spectrum  obtained  from  the  planar  detectors  included  7-rays 
with  energies  in  the  range  0  <  E1  <  1.0  MeV.  Figure  23  shows  the  summed  planar 
spectrum  over  the  energy  range  10  <  E1  <  900  keV,  outside  of  which  there  are  no 
significant  spectral  features. 


Figure  23.  Summed  planar  detector  spectrum  plotted  by  gf3,  showing  the  energies 
of  selected  prominent  peaks  from  rhenium  isotopes  and  sources  of  background.  The 
horizontal  axis  identifies  the  7-ray  energy  in  keV  according  to  the  energy  calibration 
applied.  The  vertical  axis  is  the  total  number  of  counts  per  channel,  where  the  con¬ 
version  gain  results  in  8.192  channels  per  keV.  The  134.2  keV  peak  from  187Re  and  the 
210.7  keV  peak  due  to  186Re  are  identified  to  show  the  effects  of  gating  on  the  different 
neutron  energies. 
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18'Re(n, n'7)18'Re  reaction  channel. 


The  7-rays  from  the  18'Re(n,  n/7)187Re  reaction  channel  were  selected  by  gating 
on  neutron  energies  in  the  range  of  1.0  MeV  to  10.0  MeV.  The  peaks  in  the 
resulting  spectra  are  cataloged  in  Table  2. 

Table  2.  Peaks  observed  in  planar  detectors  with  1.0  <  En  <  10.0  MeV  attributed  to 
18 'Re.  These  included  1  not  previously  reported  in  the  literature  for  18 'Re  (identified 
in  red). 


Energy  (keV) 

Fit  Error 

Database  Energy  (keV) 

Area  (counts) 

Error 

72.58 

0.01 

72.002 

33212 

590 

134.47 

0.00 

134.247 

155581 

473 

170.05 

0.00 

168.5 

72892 

376 

182.48 

0.00 

182.3 

56029 

350 

205.15 

0.01 

204.9 

24194 

293 

214.95 

0.01 

214.8 

19729 

279 

236.54 

0.01 

236.4 

15704 

266 

304.15 

0.02 

304 

7590 

237 

374.78 

0.02 

374.7 

7261 

225 

404.26 

0.03 

404 

4984 

214 

441.60 

0.05 

441.3 

2446 

444 

479.44 

0.02 

479.53 

9047 

222 

586.48 

0.04 

586.3 

6057 

549 

588.83 

0.02 

589.06 

10463 

666 

618.11 

0.02 

618.37 

9713 

250 

647.27 

0.04 

647.3 

4838 

204 

682.26 

0.05 

682.34 

3637 

224 

685.64 

0.03 

685.81 

6645 

263 

719.31 

0.07 

None  assigned 

3376 

218 

745.17 

0.09 

745.21 

2866 

199 

772.60 

0.06 

772.87 

4101 

191 

843.18 

0.04 

842.1 

11552 

518 

846.57 

0.09 

844.7 

4852 

428 
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18'Re(n,  2n7)186Re  Reaction  Channel. 


The  7-rays  from  the  18'Re(n,  2n7)186Re  reaction  channel  were  selected  by  gating 
on  neutron  energies  in  the  range  of  10.0  MeV  to  25.0  MeV.  The  peaks  in  the 
resulting  spectra  are  cataloged  in  Table  3. 

Table  3.  Peaks  observed  in  planar  detectors  with  10.0  <  En  <  25.0  MeV  attributed  to 
186Re.  These  included  4  not  previously  reported  in  the  literature  for  186Re  (identified 
in  red). 


Energy  (keV) 

Fit  Error 

Database  Energy  (keV) 

Area  (counts) 

Error 

144.37 

0.01 

143.919 

10046 

287 

210.83 

0.02 

210.685 

4146 

141 

232.41 

0.07 

232.1 

1224 

123 

266.85 

0.02 

None  assigned 

5394 

137 

316.54 

0.05 

316.473 

1577 

106 

354.31 

0.09 

None  assigned 

818 

96 

381.24 

0.06 

None  assigned 

1239 

95 

413.54 

0.06 

413.21 

527 

87 

492.64 

0.12 

None  assigned 

1424 

243 
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Coaxial  Detector  Projection  Spectra. 


The  projection  spectrum  obtained  from  the  coaxial  detectors  included  7-rays 
with  energies  in  the  range  0  <  Ey  <  4.0  MeV.  Figure  24  shows  the  summed  planar 
spectrum  over  the  energy  range  50  <  E1  <  3750  keV,  outside  of  which  there  are  no 
significant  spectral  features. 


Figure  24.  Summed  coaxial  detector  spectrum  plotted  by  gf3,  showing  the  energies 
of  selected  prominent  peaks  from  rhenium  isotopes  and  sources  of  background.  The 
horizontal  axis  identifies  the  7-ray  energy  in  keV  according  to  the  energy  calibration  ap¬ 
plied.  The  vertical  axis  is  the  total  number  of  counts  per  channel,  where  the  conversion 
gain  results  in  2.048  channels  per  keV.  No  peaks  are  identified  at  energies  greater  than 
1.1  MeV  due  to  energy  calibration  errors  that  resulted  from  ADC  gain  nonlinearities. 
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18'Re(n, n'7)18'Re  Reaction  Channel. 


The  7-ray  peaks  from  the  18'Re(n,  n77)18'Re  reaction  are  listed  in  Table  4. 

Table  4.  Peaks  observed  in  coaxial  detectors  with  1.0  <  En  <  10.0  MeV  attributed  to 
18 'Re.  These  included  4  not  previously  reported  in  the  literature  for  18 'Re  (identified 
in  red). 


Energy  (keV) 

Fit  Error 

Database  Energy  (keV) 

Area  (counts) 

Error 

134.76 

0.01 

134.247 

70187 

433 

170.20 

0.01 

168.5 

48871 

401 

182.57 

0.01 

182.3 

46003 

396 

205.19 

0.02 

204.9 

26611 

366 

214.90 

0.02 

214.8 

21064 

359 

236.43 

0.04 

236.4 

11405 

351 

304.00 

0.09 

304 

8975 

405 

374.33 

0.06 

374.7 

12350 

382 

404.07 

0.07 

404 

11353 

371 

440.78 

0.00 

441.3 

6365 

657 

454.86 

0.28 

454.92 

2483 

334 

478.95 

0.06 

479.53 

10730 

338 

492.20 

0.20 

491.2 

3299 

320 

516.57 

0.05 

518.6 

13566 

336 

585.94 

0.05 

586.3 

7199 

296 

588.23 

0.03 

589.06 

14569 

346 

617.47 

0.04 

618.37 

17871 

507 

635.96 

0.19 

635.8 

3268 

355 

646.74 

0.05 

647.3 

12222 

351 

659.84 

0.07 

None  assigned 

7375 

310 

681.76 

0.06 

682.34 

4313 

535 

685.04 

0.05 

685.81 

7773 

585 

718.97 

0.08 

None  assigned 

4335 

536 

744.45 

0.09 

745.21 

5831 

432 

771.71 

0.07 

772.87 

8456 

520 

843.46 

0.04 

842.1 

26194 

766 

960.96 

0.10 

960.17 

5133 

591 

992.45 

0.07 

None  assigned 

5964 

224 

1014.50 

0.02 

None  assigned 

25076 

275 
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18'Re(n,  2n7)186Re  Reaction  Channel. 


The  7-rays  from  the  18'Re(n,  2n7)186Re  reaction  channel  were  selected  by  gating 
on  neutron  energies  in  the  range  of  10.0  MeV  to  25.0  MeV.  The  peaks  in  the 
resulting  spectra  are  cataloged  in  Table  5. 

Table  5.  Peaks  observed  in  coaxial  detectors  with  10.0  <  En  <  25.0  MeV  attributed  to 
186Re.  These  included  3  not  previously  reported  in  the  literature  for  186Re  (identified 
in  red). 


Energy  (keV) 

Fit  Error 

Database  Energy  (keV) 

Area  (counts) 

Error 

144.58 

0.04 

143.919 

4062 

170 

210.78 

0.04 

210.685 

3612 

169 

232.32 

0.10 

232.1 

1498 

160 

251.84 

0.11 

251.841 

1396 

163 

266.58 

0.03 

None  assigned 

4511 

227 

316.49 

0.15 

316.473 

3083 

215 

380.94 

0.16 

None  assigned 

2470 

190 

390.30 

0.19 

390.91 

2032 

186 

629.65 

0.13 

None  assigned 

1258 

180 

Peaks  Above  1.1  MeV. 

Other  unassigned  peaks  with  E, y  >1.1  MeV  existed  in  the  coaxial  detector 
spectra  in  both  reaction  channels,  however,  these  peaks  were  not  attributed  to  the 
rhenium  isotopes  in  this  work.  This  was  due  to  the  fact  that  based  on  certain  7-rays 
in  the  ENSDF,  the  energy  calibration  at  7-ray  energies  above  1.1  MeV  was  observed 
to  be  in  error  greater  than  5  keV.  Apparent  nonlinearities  in  the  coaxial  detector 
gain  resulted  in  difficulty  identifying  new  transitions  from  known  sources  of 
background,  primarily  those  7-rays  due  to  (n,  xny)  reactions  in  germanium  isotopes. 
This  anomalous  behavior  persisted  despite  numerous  attempts  to  calibrate  the 
coaxial  detector  spectra  using  polynomials  of  orders  1  through  4  to  fit  the  peaks  in 
the  152Eu  source  spectrum. 
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Common  Spectral  Features. 


Both  coaxial  and  planar  detector  spectra  include  X-rays  due  to  rhenium  and 
bismuth,  the  electron-positron  annihilation  peak  at  511  keV,  and  neutron  induced 
pulses  due  to  the  direct  interaction  of  neutrons  with  the  germanium  and  bismuth 
nuclei  in  the  detectors. 


X-rays. 

X-rays  from  both  rhenium  and  bismuth  (present  in  the  BGO  suppression 

shields)  are  evident  in  both  the  coaxial  and  planar  detector  spectra.  These  X-rays 

originate  from  internal  conversion  in  the  target  or  detection  medium,  because  an 

X-ray  cascade  follows  the  emission  of  an  inner  atomic  shell  electron.  The  X-rays  are 

generally  more  intense  than  the  7-rays  due  to  neutron- induced  reactions,  so  the 

61.14  keV  Ka7  X-ray  from  rhenium  in  particular  was  used  as  the  primary 

low-energy  reference  peak  in  the  gain  matching  procedure  for  the  planar  detectors. 

Table  6.  Energies  of  rhenium  and  bismuth  X-rays  evident  in  both  planar  and  coaxial 
detector  spectra.  Values  shown  here  are  those  published  in  the  Lawrence  Berkeley  Lab 
X-ray  Data  Booklet  [42]. 


Element,  X-ray 

Energy  (keV) 

Rhenium,  Ka2 

59.7179 

Rhenium,  Koi 

61.1403 

Rhenium,  K/3| 

69.31 

Rhenium,  K/T 

71.232 

Bismuth,  K a2 

74.8148 

Bismuth,  K«! 

77.1079 

Bismuth,  K/3i 

87.343 

Bismuth,  K/32 

89.83 
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511  keV  annihilation  peak. 


The  511  keV  peak  present  in  all  the  spectra  is  clue  to  electron-positron 
annihilation  following  pair  production  in  the  detection  medium.  Gamma  rays  with 
energies  greater  than  1.022  MeV  can  interact  in  the  detector  and  result  in  the 
production  of  an  electron-positron  pair.  This  process  must  occur  within  the 
Coulomb  field  of  a  nucleus,  and  the  cross  section  for  pair  production  generally 
increases  with  Z2,  where  Z  is  the  atomic  number  of  the  detection  material  [43].  The 
electron  and  positron  depart  the  location  where  they  are  created  with  some  forward 
momentum,  dependent  on  the  energy  of  the  y-ray  from  which  they  originate.  The 
positron,  which  travels  a  short  distance  inside  the  detection  medium  before 
encountering  an  electron,  annihilates  with  the  creation  of  two  photons.  Peaks  at 
511  keV  and  1022  keV  are  evident  in  the  HPGe  spectra,  indicative  of  when  one  or 
both  of  these  two  photons  are  detected.  As  these  photons  are  often  detected  at 
distinct  times  according  to  the  time  resolution  of  the  HPGe  detector,  the  511  keV 
peak  is  significantly  more  prominent  than  the  1022  keV  sum  peak  [43]. 

Because  of  the  momenta  of  both  the  positron  and  electron  upon  annihilation, 
the  relative  velocity  of  the  particles  is  generally  some  significant  fraction  of  the 
speed  of  light  c,  and  the  annihilation  peak  is  Doppler-broadened  as  a  result. 
Compared  to  other  y-ray  peaks  in  the  ADC  spectra  from  the  planar  and  coaxial 
detectors  in  this  experiment,  the  511  keV  peak  is  significantly  wider.  Since  the 
511  keV  peak  was  clearly  evident  in  all  the  spectra,  it  was  used  for  gain  matching, 
and  had  to  be  carefully  fitted  to  account  for  its  unique  Doppler-broadened  shape. 

Neutron-induced  pulses. 

Fast  neutron  interactions  in  the  germanium  detection  medium  result  in 
characteristically  wide  ramp-shaped  peaks,  primarily  at  y-ray  energies  of 
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596.06,  609.0  and  693.4  keV.  The  first  2  of  these  3  peaks  are  due  to  (n,n' 7) 
reactions  in  '4Ge,  while  the  third  is  due  to  an  (n,  n7 7)  reaction  in  ,2Ge.  When  a  fast 
neutron  scatters  off  a  germanium  nucleus,  energy  goes  into  both  nuclear  recoil  and 
excitation  energy,  the  latter  of  which  is  released  in  7-decay.  Electron-hole  pairs  in 
the  detector  are  created  both  by  the  7-ray  and  the  recoiling  germanium  nucleus, 
which  broadens  the  peak  on  the  high-energy  side.  Other  background  contributions 
included  neutron-induced  pulses  from  209Bi  and  19F  in  the  BGO  detectors  and 
sample  holder,  respectively. 

4.4  7  —  7  Coincidence  Matrix  Analysis 

The  TScan  presort  routine  sort .  c  was  modified  to  include  definitions  for  two 
7  —  7  coincidence  matrices  gated  on  the  neutron  energies  in  the  ranges 
1.0  <  En  <  10  MeV  and  10  <  En  <  25  MeV  to  select  out  those  7-ray  events  from 
the  (n,  n'y)  and  (n,  2n 7)  reaction  channels,  respectively.  The  resulting  .  grat  matrix 
hie  was  converted  to  the  RadWare  .mat  coincidence  matrix  format  using  the 
gmt2mat  subroutine. 

RadWare  escl8r  program. 

The  escl8r  program  included  with  the  RadWare  distribution  reads  symmetric 
2-dimensional  7  —  7  coincidence  matrices  and  allows  the  user  to  gate  on  peaks  in  the 
spectra  to  see  other  7-rays  in  coincidence  with  the  selected  7-ray.  In  this  way, 
escl8r  can  be  used  to  deduce  more  information  about  the  level  scheme  for  a  given 
nucleus  than  can  be  found  from  an  analysis  of  only  singles  7-ray  spectra. 
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Search  of  Coincidence  Data  for  New  Transitions. 


The  relevant  7  —  7  coincidence  matrices  were  searched  for  evidence  of  7-ray 
cascades  involving  each  of  the  new  transitions  discovered  in  the  (n,  n' 7)  and  (n,  2nr)) 
reaction  channels  described  in  Section  4.3.  However,  at  each  new  transition  energy, 
gating  on  the  peak  resulted  in  a  complete  absence  of  other  7-rays  in  coincidence 
with  the  selected  peak  energy.  This  observation  implies  a  possibility  that  none  of 
the  new  transitions  proposed  in  this  work  are  part  of  a  7-ray  cascade,  or  that  there 
is  an  absence  of  coincidence  data  from  which  one  can  make  any  determination  about 
7-ray  cascades  in  the  186Re  and  187Re  nuclei.  An  examination  of  other  prominent 
peaks  in  the  spectra,  such  as  the  134  keV  peak  in  187Re,  supports  the  second 
hypothesis. 

4.5  New  7-ray  Transitions  in  187Re 

Tables  2  and  4  include  4  peaks  attributed  to  18'Re  not  found  in  the  ENSDF  or 
XUNDL  nuclear  databases. 

659.84(7)  keV. 

The  659.84  keV  peak  observed  in  the  coaxial  spectrum  has  a  low  intensity.  It 
was  not  discernible  above  background  in  the  planar  spectrum,  because  the  planar 
detectors  have  a  lower  detection  efficiency  than  the  coaxial  detectors  at  7-ray 
energies  above  ~  400  keV.  This  peak  is  possible  evidence  of  an  El  transition 
between  the  Jn  =  (1/2”, 3/2+, 5/2”)  level  at  1661(6)  keV  to  the  Jn  =  (5/2”,7/2+) 
level  at  1000.93(12)  keV,  which  would  have  a  calculated  energy  of  658(6)  keV  per 
the  ENSDF.  The  existence  of  the  1000.93  keV  level,  however,  is  listed  as  uncertain 
in  the  ENSDF. 
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719.1(2)  keV. 

Peaks  appear  at  energies  of  718.97(8)  keV  in  the  coaxial  spectrum  and 
719.31(7)  keV  in  the  planar  spectrum,  which  average  to  719.1(2)  keV.  These  peaks 
may  indicate  the  existence  of  an  Ml  or  E 2  transition  between  the 
Jn  =  (3/2+,5/2+)  level  at  1230.10(4)  keV  to  the  J77  =  l/2+  level  at  511.768(7)  keV. 
Such  a  transition  would  have  a  calculated  energy  of  718.35(4)  keV  per  the  ENSDF, 
implying  a  discrepancy  of  0.8  keV  between  the  observed  and  calculated  values  for 
the  energy.  The  ENSDF  shows  a  E 2  transition  from  the  511.768  level  directly  to  the 
ground  state,  but  there  is  no  such  transition  evident  in  the  7  —  7  coincidence  data. 

An  alternative  possibility  is  that  the  peak  at  719.1  keV  is  due  to  a  previously 
unobserved  transition  from  the  uncertain  718.73(4)  keV  level  (spin-parity 
unassigned)  to  the  ground  state  (J77  =  5/2+). 

992.45(7)  keV. 

A  992.45  keV  peak  was  observed  in  the  coaxial  detector  spectrum.  It  is  possible 
that  it  is  evidence  of  a  transition  from  the  J77  =  (9/2_)  level  at  1200(3)  keV  in  the 
ENSDF  to  the  Jn  =  9/2”  level  at  206.252(7)  keV,  which  has  a  calculated  energy  of 
994(3)  keV.  However,  such  a  transition  would  imply  that  the  spin-parity  assignment 
of  J77  =  (9/2_)  for  the  upper  level  is  in  error,  as  monopolc  ( L  =  0)  transitions  in 
which  a  single  photon  is  emitted  are  forbidden  by  parity  selection  rules  [14],  The 
absence  of  a  coincident  206.25  keV  7-ray  in  the  data  is  explained  by  the  fact  that 
the  206.25  keV  level  in  187Re  has  a  lifetime  of  552.3  ns,  and  decay  from  this  level  to 
the  ground  state  is  highly  converted  (a  =  3.35). 
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1014.50(2)  keV. 

The  1014.50(2)  keV  peak  observed  in  the  coaxial  spectrum  leads  to  the 
possibility  of  a  transition  between  the  level  at  1220.80(25)  keV  (spin-parity 
unassigned)  to  the  J n  =  9/2~  level  at  206.252(7)  keV.  Such  a  transition  would  have 
a  calculated  energy  of  1014.55(25)  keV.  As  in  the  case  of  the  992.52  keV  transition 
above,  no  206.252  keV  y-ray  is  evident  in  the  coincidence  data,  likely  due  to  the 
lifetime  of  the  206.25  keV  level  and  internal  conversion. 

4.6  New  y-ray  Transitions  in  186Re 

Noted  in  Tables  3  and  5  are  5  peaks  attributed  to  186Re  that  are  not  included  in 
the  ENSDF  or  XUNDL  nuclear  databases,  or  the  associated  literature.  Of  greatest 
significance  is  the  266.7  keV  y-ray  observed  in  the  planar  spectrum  that  is  expected 
to  directly  feed  the  2  x  105y  isomer.  Each  transition  is  described  in  greater  detail  in 
the  following  sections. 

The  level  schemes  of  184Re  and  186Re  show  certain  similarities  at  low  energies, 
owing  to  the  fact  that  both  isotopes  have  long-lived  isomers  in  the  140  -  190  keV 
energy  range  with  matching  fully  aligned  configurations  [1,  31].  The 

known  level  structure  above  the  188  keV  isomer  in  184Re  provided  motivation  for 
the  assignment  of  transitions  that  feed  the  isomer  in  186Re.  In  184Re,  three  primary 
transitions  that  feed  the  isomer  are  known;  these  are  shown  in  Figure  25  alongside 
known  levels  and  proposed  transitions  in  186Re. 

266.7(2)  keV. 

The  266.7(2)  keV  y-ray  is  the  average  energy  of  peaks  at  266.85(2)  keV  and 
266.58(3)  keV  in  the  planar  and  coaxial  spectra,  respectively.  This  y-ray  is  a 
candidate  for  a  transition  that  directly  feeds  the  isomer  at  149(7)  keV  from  the  level 
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Figure  25.  Extracts  from  the  184Re  and  186Re  level  schemes,  showing  similarities  in  the 
spin-parity  assignments  for  levels  in  the  two  nuclei,  and  proposed  transitions  identified 
in  this  experiment  shown  in  red.  The  629.65  keV  transition  is  dashed  to  highlight  the 
fact  that  this  is  one  of  two  possibilities  for  a  transition  of  this  energy.  Adapted  with 
permission  from  [1]. 


at  414.9(5)  keV  proposed  in  [26].  Wheldon  and  others  state  that  the  414.9  keV  level 
observed  in  the  group’s  experiment  involving  the  population  of  186Re  levels  via  the 
18'Re(p,  d)186Re  reaction  may  have  a  spin-parity  of  Jn  =  7_,  though  this  assignment 
is  only  tentative.  Here,  the  similarity  between  the  levels  in  184Re  and  186Re 
motivates  a  tentative  assignment  of  J 77  =  (9_)  for  the  spin-parity  of  the  414.9  keV 
level.  Furthermore,  existence  of  this  transition  implies  an  energy  of 
414.9(5)  —  266.7(2)  =  148.2(5)  keV  for  the  isomer,  a  significant  improvement  over 
the  149  ±  7  keV  value  originally  deduced  by  Seegmiller,  et  al.  [7] . 


354.31(9)  keV. 

The  low-intensity  354.31  keV  peak  observed  in  the  planar  spectrum  was  not 
discernible  above  background  in  the  coaxial  detector  spectrum.  This  peak  is  likely 
indicative  of  an  El  transition  between  the  J71"  =  (4)+  level  at  500.722(16)  keV  to  the 
Jn  =  (3)~  level  at  146.274(4)  keV,  which  would  have  a  calculated  energy  of 
354.448(16)  keV  per  the  ENSDF. 
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381.1(2)  keV. 

Peaks  appear  at  energies  of  381.24(6)  keV  in  the  coaxial  spectrum  and 
380.94(16)  keV  in  the  planar  spectrum,  which  average  to  381.1(2)  keV.  These  peaks 
likely  represent  evidence  of  a  El  transition  from  the  Jn  =  (5)+  level  at 
559.976(9)  keV  in  the  ENSDF  to  the  Jn  =  (6)~  level  at  180.1(5)  keV  proposed  by 
Wheldon  and  others  [26].  This  transition  would  have  a  calculated  energy  of 
379.9(5)  keV  per  the  ENSDF. 

492.64(12)  keV. 

The  492.64  keV  peak  observed  in  the  planar  spectrum  was  not  discernible  above 
background  in  the  coaxial  spectrum.  This  peak  is  likely  indicative  of  an  Ml 
transition  between  the  J77  =  (2_,3_)  level  at  1069.8  keV  to  the  J 77  =  (2“)  level  at 
577.732(16)  keV,  which  would  have  a  calculated  energy  of  492.077(16)  keV  per  the 
ENSDF.  Existence  of  such  a  transition  would  give  credence  to  a  tentative 
assignment  of  J71"  =  (3_)  to  the  1069.8  keV  level,  as  parity  selection  rules  prohibit 
monopole  ( L  =  0)  transitions  in  which  a  single  photon  is  emitted  [14], 

629.65(13)  keV. 

Motivated  by  the  184Re  level  scheme  shown  in  Figure  25,  the  629.65  keV 
transition  observed  in  the  18'Re(n,  2n7)186Re  channel  is  possibly  a  transition  from 
the  774.2(15)  keV  level  observed  by  Wheldon  and  others  that  directly  feeds  the 
2  x  105  y  isomer.  This  level  is  not  strongly  populated  in  the  187Re(p,  d)186Re 
reaction  per  Wheldon  and  others,  and  no  proposed  spin-parity  is  given  in  the 
group’s  2009  paper  [26].  Per  the  same  reasoning  as  in  the  case  of  the  267  keV  peak, 
a  spin-parity  of  J77  =  9-  is  likely  for  this  state,  as  shown  in  Figure  25.  Existence  of 
this  transition  would  imply  an  energy  of  774.2(15)  —  629.65(13)  =  144.6(15)  keV  for 
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the  isomer,  which  differs  from  the  148.2(5)  keV  value  deduced  above  by  more  than 
one  standard  deviation. 

Alternatively,  the  observed  629.65  keV  transition  could  be  evidence  of  an  Ml 
transition  from  the  =  (1“)  level  at  689.3  keV  to  the  Jn  =  (2)_  level  at 
59.010  keV,  which  would  have  a  calculated  energy  of  630.29  keV. 
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V.  Future  Work 


The  results  claimed  in  this  work  represent  the  conclusions  from  an  initial 
analysis  of  the  data.  Remaining  tasks  in  the  data  analysis  include  calculating  7-ray 
excitation  functions  and  transition  cross  sections  for  new  transitions,  and  comparing 
these  cross  sections  with  those  calculated  using  nuclear  reaction  codes. 

5.1  Excitation  Functions 

The  7-ray  excitation  function  Xy(En)  is  defined  as  the  7-ray  yield  at  a  given 
energy  divided  by  the  neutron  flux: 


v,(«  =  (14) 

Here,  J7  is  7-ray  intensity  in  photons/s  and  <f>  is  the  neutron  flux  in 
neutrons/MeV-s.  The  excitation  function  for  a  given  7-ray  is  proportional  to  the 
partial  7-ray  cross  section,  as  described  in  Section  5.2.  Obtaining  the  7-ray 
excitation  functions  for  the  new  transitions  claimed  in  this  work  will  provide  a  basis 
for  an  initial  comparison  with  theoretical  models.  These  models,  described  in 
further  detail  in  Section  5.3,  are  capable  of  producing  partial  7-ray  cross  sections,  so 
the  shapes  of  the  curves  can  be  compared.  Based  on  the  goodness  of  fit  between  the 
shapes  of  the  experimentally  obtained  7-ray  excitation  functions  and  the  partial 
7-ray  cross  sections  produced  by  the  models,  it  may  be  possible  to  deduce  transition 
multipolarities  and  spin-parity  assignments  for  certain  levels. 

5.2  Transition  Cross  Section  Calculations 

Gamma-ray  excitation  functions,  which  are  defined  as  the  7-ray  intensity  /7  as  a 
function  of  incident  neutron  energy  En  according  to  Equation  14,  will  be  used  to 
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develop  cross  sections  for  the  newly-discovered  transitions.  The  transition  cross 
section  cr(En )  for  a  particular  7-ray  energy  E1  is  related  to  the  excitation  function 
Xy(En)  through  the  following  relation  [32,  34]: 


a(En)  =  X^(En) 


DeadTime7  •  (1  +  a)  ■  C7(En) 
DeadTime$  •  e(Ey)  ■  t 


■N. 


(15) 


Here, 

•  DeadTime7;$  is  the  dead  time  correction  for  the  HPGe  and  fission  chamber 
detectors,  respectively, 

•  a  is  the  conversion  coefficient  for  the  transition, 

•  C1(En )  is  an  angular  correction  factor, 

•  e(E1)  is  the  total  efficiency  of  the  detector  array  at  the  photon  energy  E1: 

•  t  is  the  areal  density  of  the  target,  in  atoms/barn,  and 

•  N  is  a  factor  of  normalization  obtained  by  measuring  the  excitation  function 
for  a  particular  transition  in  56Fe. 

The  methods  for  obtaining  the  7  factors  listed  above  necessary  to  arrive  at  the 
transition  cross  section  from  the  excitation  function  are  explained  in  the  following 
sections. 


Dead  Time  Corrections. 

Both  the  7-ray  and  fission  chamber  ADC  spectra  must  be  corrected  for  dead 
time  to  obtain  the  number  of  counts  in  real  time  from  each  of  the  detector  types. 
This  is  critical  to  being  able  to  discern  the  actual  neutron  flux  and  7-ray  yields. 
The  process  of  correcting  for  dead  time  involves  comparing  the  detector  live  times 
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against  outputs  from  the  clock  timing  hardware,  both  of  which  are  sent  to  scaler 
modules  and  incorporated  into  the  MIDAS  data  acquisition  system.  The  scaler  data 
needed  to  calculate  the  fission  chamber  and  HPGe  detector  dead  time  are  contained 
in  the  ROOT  histograms  that  were  produced  by  the  online  MIDAS  analyzer  and 
saved  to  disk  following  the  experiment. 

Internal  Conversion  Coefficients. 

Internal  conversion  coefficients  for  some  transitions  are  published  in  the  ENSDF, 
though  a  search  of  the  database  identifies  many  transitions  for  which  no  such 
measurement  of  a  has  been  made.  This  is  particularly  true  in  the  case  of  isotopes 
that  have  not  been  the  focus  of  much  research,  such  as  186Re  [10].  When  no 
measurement  of  the  conversion  coefficient  for  a  particular  transition  exists  in  the 
database,  theoretical  models  can  be  used  to  calculate  a.  The  online  conversion 
coefficient  calculator  Brlcc,  developed  at  the  Australian  National  University,  is  one 
such  model  [23].  Given  the  atomic  number  Z,  transition  energy  E1  and  transition 
multipolarity,  the  Brlcc  code  applies  a  Dirac-Fock  computational  model  to  account 
for  electron  vacancies  in  the  process  of  internal  conversion.  It  produces  conversion 
coefficients  that  generally  match  experimental  data  to  within  several  percent  [44], 

Detector  Array  Efficiency  Calibration. 

The  HPGe  detectors  in  the  GEANIE  array  were  calibrated  individually  for 
efficiency  using  a  variety  of  calibration  sources,  a  process  that  took  place  during  the 
run  cycle. 
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Angular  Correction  for  Efficiency. 


The  absolute  7-ray  detection  efficiency  of  detectors  in  the  GEANIE  array  is 
subject  to  angular  variation  dne  to  the  inherent  multipolarity  of  the  7-rays  emitted 
from  the  target  and  effects  of  self- absorption  within  the  target.  The  product  of 
these  two  individual  factors  is  the  angular  correction  factor  C1  in  Equation  15. 

Self-absorption  effects. 

A  certain  percentage  of  7-rays  produced  in  the  target  are  absorbed  by  the  target 
before  they  can  be  detected.  Assuming  excited  nuclear  states  are  produced 
according  to  a  uniform  distribution  within  the  target,  7-rays  emitted  from  target 
nuclei  must  pass  through  varying  amounts  of  target  material  before  they  can  be 
detected.  Correcting  for  this  effect  involves  an  Monte  Carlo  Neutral  Particle 
(MCNP)  simulation  of  the  target  and  detector  array  [39,  45]. 

Multipolarity  effects. 

The  angular  distribution  of  7-rays  produced  in  a  target  relative  to  the  beam 
direction  is  described  by  the  differential  cross  section  ^|,  which  is  a  superposition  of 
even  Legendre  polynomials  Pn  (sin#),  per  the  equation 

^  (o)  =  ^  anPn  (sin0)  (16) 

n  even 

where  <j0  is  the  total  cross  section,  6  is  the  scattering  angle,  and  the  an  are  the 
Legendre  polynomial  coefficients  [34,  39].  The  total  cross  section  is  calculated  from 
the  total  detector  array  efficiency,  and  the  individual  detector  locations  listed  in 
Table  1  provide  the  scattering  angle  9.  A  code  known  as  AVALANCHE  is  available 
for  performing  these  angular  correction  calculations  given  a  set  of  parameters 
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defining  the  detector  geometry. 


Normalization  Using  56Fe. 

The  normalization  constant  N  in  Equation  15  relates  the  excitation  function  to 
the  transition  cross  section  a.  This  normalization  constant  is  obtained  by  measuring 
the  7-ray  yield  for  the  846.8  keV  transition  in  56Fe  from  the  reaction 
56Fe(n,  n/7)56Fe.  The  cross  section  for  this  reaction  is  well  known  and  its  value  at 
the  neutron  energy  En  =  14.5  MeV  is  used  as  the  normalization  constant. 
Application  of  this  normalization  factor  is  necessary  to  correct  for  the  differences  in 
the  sizes  of  the  beam  spot  and  the  target. 

5.3  Reaction  Modeling 

Two  primary  computational  codes  are  available  that  model  partial  7-ray  cross 
sections,  equal  to  the  transition  cross  section  divided  by  (1  +  a),  where  a  is  the 
internal  conversion  coefficient  for  the  transition  in  question.  These  codes,  COH3  and 
TALYS  1.6,  are  commonly  used  to  validate  experimental  data  obtained  at  GEANIE 
[34,  35].  COH3,  which  was  developed  at  LANSCE,  calculates  partial  7-ray  cross 
sections  for  neutron- induced  reactions.  Both  codes  have  been  used  to  validate 
transition  cross  sections  obtained  from  recent  GEANIE  experiments  [34,  35,  37]. 

COH3. 

COH3  is  a  C++  code  developed  at  LANSCE  by  T.  Kawano  for  the  purpose  of 
modeling  nuclear  reaction  cross  sections  in  mid-  to  liigh-Z  elements.  It  supersedes 
the  earlier  GNASH  code,  built  for  the  same  purpose  and  used  to  validate 
experimental  cross  sections  prior  to  2008.  COH3  is  capable  of  producing  differential 
7-ray  cross  sections  given  neutron  energies  in  the  keV  to  MeV  energy  range,  which 
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permit  it  to  be  used  to  validate  cross  sections  obtained  in  this  experiment.  The 
theoretical  models  used  to  determine  these  cross  sections  are  optical  and 
Hauser- Feshbach  statistical  models,  and  the  code  applies  these  using  a  Monte  Carlo 
algorithm  [46,  47]. 

TALYS. 

The  TALYS  program  is  designed  to  model  a  variety  of  nuclear  reactions  in  mid- 
to  liigh-Z  elements  in  the  energy  range  1  keV  -  200  MeV.  TALYS,  which  is 
open-source  software,  was  developed  as  part  of  a  collaboration  between  the  French 
Alternative  Energies  and  Atomic  Energy  Commission  (CEA)  and  the  Nuclear 
Research  and  Consultancy  Group  (NRG)  in  the  Netherlands.  The  code  is  used  to 
simulate  nuclear  reactions  for  the  purposes  of  analyzing  experiments  and  generating 
nuclear  data,  including  cross  sections,  in  applications  where  there  is  insufficient 
experimental  data  [48].  The  code  relies  on  a  number  of  nuclear  reaction  models, 
including  the  Hauser-Feshbach  model  that  is  also  used  in  the  COH3  nuclear  reaction 
code.  TALYS  is  a  FORTRAN-based  program  that  generates  cross  sections  and 
reaction  Q-values  given  user  input  consisting  of  a  deck  of  reaction  parameters. 

5.4  Examination  of  Data  from  natHe  Runs 

Similarities  between  the  level  structures  of  184Re  and  186Re  suggest  it  is  possible 
that  the  excitation  functions  for  the  transitions  that  feed  the  J71"  =  8+  isomer  in 
184Re  may  be  similar  to  those  that  feed  the  isomer  in  186Re.  To  test  this  hypothesis, 
24  hours  of  irradiation  of  sample  of  natural  rhenium,  which  contains  37.4%  185Re, 
was  accomplished  on  8  -  9  December.  The  data  obtained  from  these  experimental 
runs  will  be  used  to  obtain  excitation  functions  for  y-rays  produced  by  the 
185Re(n,  n'7)185Re  and  185Re(n,  2ny)184Re  reactions. 
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Excitation  Functions  in  Thallium  Isotopes. 


An  analysis  cross  sections  for  the  new  transitions  discovered  in  186Re  will  follow 
the  methodology  used  in  a  series  of  GEANIE  experiments  in  2007  -  2008,  in  which 
previously  unknown  transitions  populating  the  J71"  =  7+  isomer  in  204T1  were 
identified.  These  experiments,  which  were  designed  to  study  the  nuclear  structure  of 
204T1,  identified  several  new  transitions  in  201T1  from  the  205Tl(n,  2n7)204Tl  reaction. 

The  existence  of  the  7+  isomer  in  204T1  was  previously  known,  although  the 
structure  above  the  isomer  was  not  [49].  Similarities  between  the  level  schemes  of 
202T1  and  204T1  led  to  the  hypothesis  that  the  structures  above  the  isomer  were 
similar  [37].  Both  202T1  and  201T1  have  a  J77  =  7+  isomer  with  similar  energies  and 
half-lives.  These  similarities  between  the  two  isotopes  suggested  that  the  levels 
above  the  isomer  in  204T1  were  similar  to  those  in  202T1. 

Excitation  functions  for  new  transitions  observed  in  the  205Tl(n,  2rr)/)204Tl 
reaction  channel  were  compared  with  those  for  the  analogous  transitions  from 
spectra  taken  from  the  203Tl(n,  2n7)202Tl  reaction.  The  results  showed  similarities 
between  the  two  sets  of  excitation  functions.  Evidence  from  modeling  and  other 
experiments  led  to  the  conclusion  that  for  transitions  with  identical  final  level  spin 
parity  assignments,  similar  excitation  function  shapes  implied  the  same  initial  level 
spin  parity.  This  led  to  the  tentative  assignment  of  Jn  =  7+,  for  the  1454  keV  level 
in  204T1  that  populates  the  8+  isomer  via  a  349.9  keV  direct  transition. 

Comparison  of  Rhenium  Excitation  Functions. 

Proceeding  as  in  the  case  of  the  thallium  isotopes,  the  data  from  the  natRe  runs 
will  be  used  to  obtain  7-ray  excitation  functions  from  the  185Re(n,  2717) 184 Re 
reaction  channel.  The  excitation  functions  for  the  known  7-ray  transitions  feeding 
the  Jn  =  8+  isomer  in  184Re  will  be  used  with  the  excitation  functions  obtained 
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from  the  187Re(n,  277-7) 186 Re  reaction  channel  to  motivate  spin-parity  assignments 
for  levels  above  the  isomer  in  186Re  [31]. 
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VI.  Conclusion 


The  significant  outcome  of  this  experiment  was  enriching  the  level  schemes  of 
18 'Re  and  186Re.  In  particular,  the  structure  above  the  2  x  105  y  isomer  in  186Re  has 
not  been  proposed  before,  and  new  7-ray  transitions  observed  in  the 
187Re(n,  2717) 186 Re  reaction  channel  may  be  evidence  of  pathways  by  which  the 
isomer  is  populated.  One  of  these  proposed  transitions  leads  to  an  improved 
estimate  of  the  excitation  energy  of  the  186mRe  isomer  with  a  significantly  smaller 
uncertainty  than  the  value  adopted  in  the  literature. 

6.1  Contributions  to  the  Level  Schemes  of  187Re  and  186Re 

The  new  transitions  proposed  here  as  a  result  of  the  preliminary  analysis  of  the 
data  from  the  November-December  2014  GEANIE  experiment  represent  a 
contribution  to  the  enrichment  of  the  level  schemes  of  both  18 'Re  and  186Re. 

Though  the  primary  focus  of  this  research  was  on  the  nuclear  structure  of  186Re, 

4  new  levels  were  also  observed  in  18 'Re.  Among  the  5  new  transitions  discovered  in 
186 Re,  one  was  identified  as  feeding  the  2  x  105  y  isomer  in  186Re.  If  validated  by 
nuclear  reaction  models,  this  will  represent  a  substantial  development  towards 
understanding  the  structure  above  this  especially  long-lived  isomer. 

The  absence  of  an  observed  cascade  to  the  ground  state  from  one  of  the  possible 
high-spin  states  above  the  isomer  would  imply  that  any  such  decay  mode  has  a  low 
intensity  relative  to  the  transitions  that  populate  the  isomer.  As  a  result,  it  is 
unlikely  an  intermediate  state  was  discovered  by  which  the  isomer  could  be 
effectively  depleted  in  the  scheme  described  in  Section  2.5. 
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6.2  Improved  Estimate  of  186mRe  Energy 


The  ENSDF  database  states  the  energy  of  the  2  x  105  y  isomer  in  186Re  is 
149  ±  7  keV  above  the  ground  state,  the  original  energy  proposed  by  Seegmiller  and 
others  in  1972  [7].  The  observed  7-ray  peak  at  267(1)  keV  in  the  187Re(n,  2n7)186Re 
reaction  channel  is  proposed  in  this  work  to  be  a  transition  to  the  isomer  state  from 
the  414.9(5)  keV  level.  This  places  the  isomer  energy  at 

414.9(5)  —  266.7(2)  =  148.2(5)  keV,  which  represents  significant  improvement  over 
the  149(7)  keV  value  in  the  evaluated  database. 

6.3  Recommendations  for  Future  Research 

Calculating  7-ray  excitation  functions  for  the  new  transitions  proposed  in  this 
work  will  permit  an  initial  comparison  of  the  experimental  data  with  transition 
cross  sections  produced  by  the  nuclear  reaction  codes  Talys  1.6  and  COH3.  After 
excitation  functions  are  obtained,  transition  cross  sections  can  be  calculated  by 
applying  corrections  for  detector  dead  time  and  efficiency,  and  these  cross  sections 
can  be  directly  validated  against  theoretical  models. 

Data  obtained  from  the  experiment  includes  the  results  of  24  hours  of  irradiation 
of  a  natural  rhenium  target.  Due  to  the  presence  of  185Re  in  the  natural  rhenium 
sample,  these  spectra  include  7-rays  produced  via  the  185Re(n,  n'7)185Re  and 
185Re(n,  2ny)184Re  reactions.  From  this  data,  7-ray  excitation  functions  from 
transitions  in  184Re  that  are  analogous  to  the  ones  proposed  for  186 Re  may  provide  a 
basis  for  proposing  spin-parity  assignments  in  186Re. 

Following  submission  of  this  thesis  document,  the  Ph.D.  dissertation  will 
continue  to  examine  the  data  obtained  from  the  GEANIE  experiment,  and  the 
dissertation  will  include  this  additional  data  analysis,  as  described  in  Chapter  V. 
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Appendix  A.  Photoexcitation  of  115mIn  Using  a  Dynamitron 

Accelerator 


Though  not  directly  related  to  the  rhenium  experiment  described  in  the  body  of 
the  thesis,  this  appendix  describes  a  demonstration  of  photoexcitation  of  nuclear 
isomers,  which  is  analogous  to  the  photodepletion  scheme  by  which  the  186mRe 
isomer  could  be  induced  to  release  its  excess  energy  once  a  suitable  intermediate 
state  is  discovered.  The  experiments  described  in  this  appendix  also  served  the 
purpose  of  checking  the  energy  calibrations  of  the  Dynamitron  accelerator  at  Wright 
State  University  (WSU)  and  the  Varian  L200  LINAC  at  the  Army  Research 
Laboratory  (ARL).  Both  of  these  accelerators  may  be  used  in  future  studies 
involving  nuclear  isomers,  so  it  is  important  that  they  are  accurately  calibrated. 

Photoexcitation  of  nuclear  isomers  by  (7,  7’)  reactions  is  possible  by  irradiation 
with  bremsstrahlung  produced  with  an  electron  LINAC,  according  to  the  theory 
outlined  in  Section  2.5  of  the  thesis.  In  the  experiment  described  in  this  appendix, 
115In  nuclei  were  excited  to  the  336  keV  isomer  115mIn  via  the  115In(7, 7')115mIn 
reaction  using  the  WSLI  Dynamitron  operated  at  1.2  MV.  115mIn  7-decays  directly 
to  the  ground  state  with  a  half  life  of  4.486  hours,  which  allowed  for  measurement  of 
the  sample  activity  after  irradiation  in  a  HPGe  detector  at  AFIT. 

The  results  of  the  WSU/AFIT  experiment  were  compared  with  those  from  a 
similar  experiment  involving  a  Varian  L200A  LINAC  operated  at  a  nominal  voltage 
of  1.0  MV  at  the  Army  Research  Laboratory  (ARL)  in  Adelphi,  MD.  A  significant 
difference  in  the  115mIn  yields  achieved  in  the  two  experiments  leads  to  the 
hypothesis  that  one  of  the  two  accelerators  operated  at  an  energy  different  from 
that  stated. 

By  modeling  the  bremsstrahlung  spectra  from  each  experiment  configuration 
using  GEANT4,  theoretical  yields  for  the  production  of  the  115mIn  isomer  were 
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calculated  from  known  cross  sections  for  the  115In(7,  7’)115mIn  reaction.  By 
comparing  the  observed  115mIn  activities  against  the  calculated  values,  it  was 
possible  to  check  the  energy  calibrations  of  the  Dynamitron  accelerator  and  the 
Varian  LINAC. 

1.1  Introduction 

The  Dynamitron  at  WSU  can  be  used  to  produce  bremsstrahlung  photons  via  a 
tungsten  converter  with  endpoint  energies  np  to  1.2  MeV.  The  advantage  of  the 
Dynamitron  is  the  relatively  high  current  (<  10  /j>A)  that  it  can  achieve  compared 
to  that  of  a  van  de  Graaf  accelerator,  as  well  as  its  continuously  tunable  energy. 

The  facility  at  WSU  is  thus  capable  of  investigating  the  population  and  depletion  of 
nuclear  isomers  by  low  energy  (7,7’)  reactions. 

Depending  on  the  particular  nucleus,  a  given  isomer  might  be  fed  by  numerous 
intermediate  states,  and  there  may  be  a  host  of  photon  energies  that  can  be  used  in 
photoexcitation  experiments.  In  the  case  of  115In,  the  first  three  IS’s  that  feed  the 
336  keV  isomer  are  the  Jn  =  7/2+  state  at  934  keV,  the  Jn  —  5/2+  state  at  941  keV 
and  the  Jn  =  ll/2+  state  at  1078  keV.  Between  1.1  MeV  and  3.0  MeV,  there  are  at 
least  five  others  [17]. 

Given  the  partial  level  scheme  shown  above,  it  is  possible  to  deduce  the 
integrated  cross  sections  for  inelastic  photon  excitation  using  Equation  8.  To  do  so, 
however,  the  lifetime  of  the  597  keV  level  is  required,  and  the  ENSDF  provides  only 
an  upper  limit  of  r  <  0.25  ns  [10].  The  natural  decay  width  T  for  this  state,  which 
is  a  function  of  the  level  lifetime  r  according  to  T  =  K/t,  is  thus  subject  to  a  large 
degree  of  variation.  As  a  result,  integrated  cross  sections  for  the 
934,  941  and  1078  keV  levels  cannot  be  accurately  calculated.  Instead,  integrated 
cross  sections  can  be  directly  measured,  the  results  of  which  have  been  published  by 
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Figure  26.  Simplified  low-energy  level  scheme  of  115In.  The  336  keV  isomer  is  fed 
by  levels  at  934  keV,  941  keV  and  1078  keV  in  the  115In(7,  7')115mIn  reaction.  The 
pathways  involving  the  934  keV  level  by  which  the  isomer  is  fed  are  highlighted  in  red. 
Figure  copyright  2001  Elsevier  B.V.  Reproduced  with  permission  from  [17]. 


Belie  and  others  for  the  photoexcitation  of  the  115mIn  isomer  [17].  The  lowest  energy 
integrated  cross  sections  for  this  reaction  are  listed  in  Table  7. 

Table  7.  The  lowest-energy  intermediate  states  (IS)  for  photoexcitation  of  the  115mIn 
isomer  and  their  associated  integrated  cross  sections  (ICS)  [17]. 


Intermediate  state  cross  sections 

Assigned  levels 

EIS  (keV) 

ICS  (ev-barn) 

Ex  (keV) 

r 

940 

0.0084  ±  (0.0004)stet.  ±  (0.0023)^. 

934 

941 

7/2+ 

5/2+ 

1085 

0.140  ±  (0.001W  ±  (0.04)^. 

1078 

11/2+ 

1449 

9/2+ 

1463 

7/2+ 

1490 

0.78  ±  (0.06)stat.  ±  (0.21)iyit. 

1487 

9/2+ 

1497 

7/2+ 
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It  is  evident  from  this  table  that  the  LIN  AC  must  be  capable  of  operating  at  a 
voltage  greater  than  940  keV  to  achieve  any  yield  of  115mIn  from  photoexcitation. 
Furthermore,  the  1085  keV  cross  section  is  greater  than  that  at  940  keV  by  a  factor 
of  17,  so  a  plot  of  yield  versus  beam  energy  would  show  a  marked  nonlinearity  at 
this  energy. 

It  is  possible  to  determine  the  yield  of  the  115mIn  isomer  created  by 
photoexcitation  from  measurements  of  the  decay  activity  of  the  indium  sample 
following  irradiation.  115In  is  nearly  stable  {t\/2  =  4.41  x  1014  y),  while  its  first 
isomeric  state  115mIn  decays  via  emission  of  a  336  keV  y-ray  to  the  ground  state 
with  a  branching  fraction  of  45.8%.  The  isomer  has  a  half  life  of  t\i 2  =  4.486  h,  so 
the  activity  is  observable  well  after  irradiation  has  ceased. 

1.2  Hypothesis 

The  Dynamitron  accelerator  at  WSU  was  used  twice  in  the  spring  of  2014  to 
irradiate  a  sample  of  natural  indium  in  separate  attempts  at  achieving  measurable 
amounts  of  115mIn.  The  y-ray  activities  were  measured  with  a  3”x3”  Nal(Tl) 
detector  following  irradiation  times  of  4.5  and  9.0  hours,  respectively,  with  an 
operating  voltage  of  1.15  MV  and  beam  current  of  5  pA.  Neither  experiment 
resulted  in  an  identifiable  336  keV  peak  in  the  measured  spectra. 

A  similar  experiment  was  conducted  at  the  Army  Research  Laboratory  in 
Adelphi,  MD,  on  5  May,  2014,  using  a  Varian  L200A  LINAC  set  at  a  nominal 
voltage  1.0  MV  and  100  /iA  beam  current.  In  this  experiment,  a  natural  indium 
sample  was  irradiated  for  a  period  of  2.0  hours  and  y-ray  activity  was  measured 
using  a  high-purity  germanium  (HPGe)  detector.  The  y-ray  spectra  from  the 
irradiated  indium  sample  in  the  ARL  experiment  clearly  showed  a  peak  at  336  keV. 

The  results  of  these  experiments  suggest  the  energy  calibration  of  either  the 


83 


Dynamitron  or  Varian  LINAC  is  in  error.  This  hypothesis  is  partially  supported  by 
the  results  of  a  calibration  of  the  Varian  LINAC  at  ARL  performed  by  Feroli  and 
others  in  2009  using  a  Compton-Hall  voltmeter.  The  researchers  concluded  from 
this  measurement  that  the  LINAC  operates  at  a  voltage  in  excess  of  1.2  MV  at  the 
nominal  1.0  MV  setting  [50].  Without  measuring  the  electron  beam  energy  or 
bremsstrahlung  photon  energy  spectra  directly,  measurements  of  photoexcitation  of 
the  115mIn  isomer  allow  one  to  verify  the  energy  calibrations  of  the  two  accelerators. 
By  comparing  observed  activities  due  to  115mIn  with  the  results  of  theoretical  yield 
calculations,  the  calibration  error  on  the  affected  accelerator  can  be  estimated.  This 
method  was  motivated  by  results  published  by  L.  Earwaker  and  D.  Weaver,  who 
used  resonant  nuclear  reactions  to  diagnose  faults  in  a  Dynamitron  at  the 
Birmingham  Radiation  Centre,  U.K.  [51]. 

1.3  Experiment 

The  two  experiments  relevant  to  this  discussion  took  place  at  ARL  on 
May  5,  2014,  and  at  WSU/AF1T  on  August  7,  2014. 

WSU/AFIT  Experiment. 

In  order  to  test  the  hypothesis,  a  measurable  yield  of  115mIn  had  to  be  achieved 
with  the  Dynamitron,  so  the  photoexcitation  experiment  was  repeated  at  WSU. 
Following  irradiation  for  3.9  hours  at  a  voltage  of  1.2  MeV  and  a  beam  current  of 
1.5  /xA,  the  sample  was  transported  to  AF1T  for  measurement  in  a  Canberra  HPGe 
detector.  This  time,  the  y-ray  spectra  taken  after  irradiation  included  a  clear  peak 
at  336  keV,  likely  due  to  the  improved  background  suppression  achieved  with  a 
copper-lined  lead  cave  surrounding  the  detector.  The  y-ray  activity  was  sufficiently 
low  (on  the  order  of  1  Bq)  that  the  peak  would  have  been  below  the  background 
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threshold  if  the  sample  were  measured  using  an  Nal(Tl)  detector  as  before. 


Dynamitron. 

The  1.5  MV  Dynamitron  at  WSU  cannot  operate  at  voltages  greater  than 
1.2  MeV,  due  primarily  to  vacuum  system  constraints.  At  1.2  MV,  the  accelerator  is 
capable  of  beam  currents  of  approximately  1.5  /J.A,  which  was  the  current  used  for 
the  indium  photoexcitation  experiment. 

Converter  and  Target. 

In  order  to  maximize  the  flux  of  bremsstrahlung  photons  incident  on  the  indium 
target,  the  converter  assembly  and  5.0  mm-thick  natural  indium  (95.7%  115In)  foil 
sample  were  positioned  at  90°  to  the  incident  electron  beam.  A  typical  Dynamitron 
target  configuration  is  shown  in  Figure  27,  with  the  converter  and  indium  foil 
mounted  against  the  copper  cold  head  that  forms  the  beam  tube  endcap.  To 
minimize  self-absorption,  a  converter  was  selected  that  consisted  of  the  thinnest 
piece  of  tungsten  (0.25  mm)  immediately  available. 


Figure  27.  Photograph  of  indium  target  and  tungsten  converter  mounted  on  the  copper 
cold  head  in  the  evacuated  Dynamitron  target  chamber. 
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To  eliminate  the  possibility  of  non-resonant  nuclear  excitation  resulting  from 
electron  scattering  in  the  indium,  an  aluminum  spacer  was  needed  between  the 
tungsten  and  indium  foils  to  absorb  electrons  that  passed  through  the  tungsten  [19]. 
The  minimum  thickness  of  aluminum  needed  was  determined  using  CASINO 
(version  3.2),  a  Monte-Carlo  simulation  program  specifically  designed  to  model  the 
trajectories  of  electrons  through  solid  materials  [52],  The  results  of  a  CASINO 
simulation  of  1000  electrons  at  1.2  MeV  (Figure  28)  were  used  to  determine  that  a 
1.8  mm-thick  piece  of  aluminum  was  sufficient  to  stop  electrons  from  scattering  into 
the  indium  target. 


Figure  28.  CASINO  simulation  to  determine  thickness  of  the  tungsten  converter  and 
aluminum  spacer  necessary  to  stop  electrons.  The  electron  beam  is  incident  from  the 
top  of  the  figure,  while  the  vertical  axis  is  penetration  depth  into  the  material.  The 
blue  lines  represent  the  paths  of  the  simulated  electrons  in  the  material,  while  the  red 
lines  are  the  characteristic  X-rays  from  tungsten. 


The  aluminum  separator  ultimately  used  in  the  converter  assembly  was  3.1  mm 
thick,  chosen  based  on  the  sizes  of  pieces  immediately  available  for  the  experiment. 


Irradiation. 


The  indium  sample  was  irradiated  for  a  combined  period  of  3  hours,  55  minutes, 
using  a  voltage  of  1.2  MeV  and  a  beam  current  of  1.5  /xA.  When  operating  at 
voltages  in  excess  of  1.0  MeV,  the  corona  current  in  the  Dynamitron  pressure  vessel 


underwent  excursions  that  required  the  beam  current  to  be  decreased  for  short 
periods  of  time.  Also,  because  of  an  inoperative  cooling  system,  the  Dynamitron 
had  be  shut  down  approximately  once  per  hour  to  cool  the  cold  head  assembly  with 
ice.  These  two  factors  prevented  the  beam  current  from  being  maintained 
continuously  at  1.5  /i A.  However,  an  average  beam-on  electron  flux  of 
2.05  x  1012  electrons/cm2 s  was  calculated  from  the  total  charge  shown  on  the 
current  integrator  at  the  cessation  of  irradiation.  The  beam-on  and  beam-off  times 
were  recorded  and  included  in  the  yield  calculations. 

Sample  Activity  Measurement. 

Following  irradiation,  the  indium  sample  was  transported  to  AFIT  for 
measurement  in  a  Canberra  HPGe  detector.  The  time  between  beam  off  in  the 
Dynamitron  and  the  beginning  of  sample  activity  measurement  was  23  minutes, 
which  was  incorporated  into  the  yield  calculations. 

The  HPGe  detectors  at  AFIT  sit  in  copper-lined  lead  caves  for  background 
suppression,  and  the  detectors  are  integrated  with  Canberra’s  proprietary  DSA-1000 
high-voltage  power  supply  /  spectroscopy  amplifier  /  multichannel  analyzer  (MCA). 
The  DSA-1000  hardware  interfaces  with  Genie  2000  software  that  runs  on  a  desktop 
PC  in  the  laboratory.  The  DSA-1000  incorporates  pulse  pile-up  rejection  (PUR) 
circuitry  that  discards  signals  due  to  pulse  pile-up  [43,  53].  Since  the  count  rates 
were  very  low  (on  the  order  of  1  s'1),  the  PUR  parameters  were  set  to  the  minimum 
allowable  values  to  avoid  pulse  rejection.  Batch  mode  processing  of  the  measured 
y-ray  spectra  was  not  necessary  for  signal  smoothing  or  further  pile-up  rejection,  as 
the  acquired  spectra  had  sufficient  definition  to  easily  discern  the  336  keV  peak  due 
to  the  presence  of  115mIn  in  the  sample. 

Before  measuring  the  indium  sample  activity,  the  HPGe  detector  was  calibrated 
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for  energy  and  efficiency  by  measuring  the  7-ray  spectrum  from  a  multinuclide 
planchette  calibration  source  for  600  seconds.  Present-day  activities  for  each 
radionuclide  in  the  source  were  calculated  and  input  in  Genic  2000  to  perform  a 
manual  efficiency  calibration,  the  results  of  which  are  shown  in  Figure  29. 


Figure  29.  Absolute  detector  efficiency  curve  resulting  from  calibration  of  the  Canberra 
HPGe  detector  with  a  multinuclide  planchette  source. 

A  least-squares  fit  was  overlaid  on  the  efficiency  data  to  determine  the  HPGe 
detector  efficiency  at  336  keV,  which  was  necessary  for  comparing  experimental  and 
theoretical  yields. 

Sample  activities  were  on  the  order  of  picocuries,  so  hourly  (3600-second  live 
time)  spectra  were  measured  to  ensure  acceptable  counting  statistics.  These  hourly 
measurements  were  obtained  over  a  period  of  7  hours  following  irradiation,  and  a 
final  spectrum  was  measured  the  following  morning,  17  hours  after  irradiation 
ceased.  The  resulting  peak  counts  were  used  to  measure  the  half-life  and  initial 
activity  of  the  1 15mIn  in  the  sample  (Figure  30). 

Agreement  between  the  measured  half  life  and  the  t\/ 2  =  4.486  hour  value  from 
literature  confirmed  that  the  decay  of  115mIn  was  the  source  of  the  336  keV  7-ray 
peak. 
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Time  after  sample  arrived  at  detector  [hours] 


Figure  30.  Hourly  counts  under  336  keV  peak  as  measured  in  HPGe  detector  for  the 
WSU/AFIT  experiment.  Fitting  an  exponential  decay  curve  to  the  hourly  peak  counts 
allowed  calculation  of  the  initial  activity  Ao  and  half  life  t1/2- 

ARL  Experiment. 

The  2  MeV  Varian  L200A  electron  LINAC  at  ARL  was  originally  designed  for 
industrial  radiography  applications.  For  this  reason,  a  tungsten  converter  was 
incorporated  in  the  LINAC  at  the  electron  beam  terminus,  similar  to  the 
configuration  used  in  the  Dynamitron.  The  higher  electron  energies  achievable  with 
the  Varian  L200A  mandate  the  incorporation  of  a  0.84  mm  thick  tungsten 
converter  [50]. 

The  target,  which  consisted  of  a  single  0.15  mm  thick  n“Hn  foil,  was  irradiated 
for  a  period  of  2.0  hours  at  an  operating  voltage  of  1.0  MV  and  current  of  100  fiA, 
followed  immediately  by  measurement  of  the  sample  activity  in  a 
background-suppressed  HPGe  detector  for  a  live  time  of  3600  seconds.  The  absolute 
detector  efficiency  at  336  keV  was  3.7% 


Comparison. 


The  3600-second  7-ray  spectra  obtained  from  the  first  hour  after  irradiation  in 
both  experiments  are  compared  in  Figure  31. 
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(a)  WSU/AFIT  Spectrum  (b)  ARL  Spectrum 

Figure  31.  First  hour  7-ray  spectra  obtained  from  measurement  of  irradiated  indium 
samples  from  (a)  WSU/AFIT  experiment  and  (b)  ARL  experiment,  both  taken  using 
background-suppressed  HPGe  detectors.  Note  the  difference  in  vertical  scales  be¬ 
tween  the  two  plots.  The  511  keV  and  1461  keV  peaks  in  (a)  are  background  due  to 
positron  annihilation  following  pair  production  and  environmental  40K.  The  small  peak 
at  392  keV  in  (b)  is  due  to  the  7-decay  of  the  t1/2  —  99.5  hour  isomer  113mIn. 


Notably,  the  ARL  spectrum  shows  53%  more  counts  under  the  336  keV  peak, 
despite  the  fact  that  the  indium  target  used  in  the  WSU/AFIT  experiment  was 
33  times  thicker  than  that  used  in  the  ARL  experiment. 


1.4  Modeling  Approach 

In  order  to  calculate  a  theoretical  yield  for  the  production  of  the  isomer  115mIn, 
one  needs  the  spectral  flux  density  ^  of  the  photon  source,  from  E  =  0  to  the 
accelerator  endpoint  energy  E0}  per  Equation  7.  The  relative  intensity  at  the 
resonant  photon  energy  for  the  reaction  is  proportional  to  the  total  flux  (a  function 
of  the  beam  current)  and  the  integrated  cross  section  of  the  appropriate  resonant 
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transition,  further  explained  in  Section  2.5  of  the  thesis.  The  spectral  flux  density 
was  simulated  using  GEANT4  instead  of  directly  measured.  Simulated  photon 
energy  spectra  for  the  WSU/AF1T  and  ARL  experimental  configurations  were  used 
to  calculate  theoretical  yields  for  photoexcitation  of  the  115mIn  isomer.  These 
calculated  yields  were  then  compared  with  those  obtained  experimentally  from 
measuring  7-ray  activities  with  HPGe  detectors. 

GEANT4. 

GEANT4  was  used  to  model  both  experiments  to  obtain  spectral  flux  densities 
at  the  indium  target  locations.  Both  GEANT4  simulations  used  the  same  step, 
track,  run  and  event  action  codes.  For  each  bremsstrahlung  photon  that  comes  from 
outside  the  indium  target  volume  and  interacts  within  the  target  volume,  these 
codes  wrote  the  photon  energy  to  a  text  hie  that  was  subsequently  used  to  generate 
an  energy  histogram  in  Matlab.  The  detector  construction  and  primary  generator 
action  codes  were  unique  to  each  experiment,  and  varied  according  to  beam  energy 
and  geometry. 

WSU/AFIT  Experiment. 

Using  the  accelerator,  converter  and  target  geometry  described  in  Section  1.3,  a 
C++  subroutine  DetectorConstruction.cc  was  written  to  model  the  WSU/AFIT 
experiment  geometry  in  GEANT4.  With  the  subroutine 

PrimaryGeneratorAction.ee  subroutine  designed  to  simulate  a  7.9  mm  diameter 
beam,  the  code  was  run  for  approximately  6  days  to  simulate  1  x  109  electrons  at 
1.2  MeV. 
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ARL  Experiment. 


The  Varian  L200A  geometry  was  modeled  in  GEANT4  to  determine  the  photon 
flux  incident  on  the  indium  foil  sample.  The  ARL  simulation  was  run  for 
approximately  7  days  to  simulate  1  x  109  electrons.  The  beam  energy  used  in  the 
simulation  was  1.3  MeV,  due  to  the  conclusion  of  Feroli  and  others  that  the  LINAC 
operating  voltage  was  “greater  than  1.2  MV”  at  the  nominal  1.0  MV  setting  [50]. 

Spectral  Flux  Densities. 

The  spectral  flux  densities  obtained  from  GEANT4  were  binned  in  1  keV  energy 
bins,  the  results  of  which  are  shown  in  the  histogram  plots  of  Figure  32.  The 
majority  of  each  plot  was  well-fitted  by  an  exponential  curve  f(x)  =  a  ■  ebx,  but  the 
tails  of  the  two  plots  (energies  >  0.9  x  E0)  were  nearly  linear.  The  tails,  which 
included  the  energies  of  the  intermediate  states,  were  therefore  fitted  using  a  linear 
function  f(x )  =  a  +  bx,  with  R2  values  of  0.9948  and  0.9932  for  the  WSU/AF1T  and 
ARL  spectra,  respectively. 

These  curves  were  used  to  determine  the  flux  at  940  keV  and  1085  keV,  which 
were  then  applied  in  the  calculation  of  isomer  yield. 

1.5  Results  and  Analysis 

The  photon  flux  at  940  keV  and  1085  keV  in  the  plots  shown  in  Figure  32  were 
used  with  published  integrated  cross  sections  from  Table  7  to  calculate  theoretical 
7-ray  activities  clue  to  the  decay  of  115mIn.  The  initial  activity  A0  in  decays/s,  or 
Bequerels  (Bq),  when  each  indium  sample  arrived  at  the  detector  following 
irradiation  is  proportional  to  the  yield  achieved  by  photoexcitation,  so  it  can  be 
used  as  a  basis  of  comparison  between  the  theoretical  and  experimental  results. 
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(a)  Dynamitron  (b)  Varian  L200A 

Figure  32.  Modeled  spectral  flux  densities  from  (a)  WSU/AFIT  and  (b)  ARL  exper¬ 
iments.  The  peaks  at  59  keV  and  67  keV  are  characteristic  X-rays  from  tungsten. 
The  figures  also  identify  the  flux  values  at  the  energies  for  which  the  integrated  cross 
sections  are  published. 


WSU/AFIT  Experiment. 


The  theoretical  and  experimental  results  from  the  WSU/AFIT  experiment 
showed  agreement  to  within  one  standard  deviation: 


Theoretical:  A0  =  0.71  ±  0.18  Bq 

Experimental:  A0  =  0.91  ±  0.12  Bq 

The  relatively  close  agreement  between  the  theoretical  and  experimental  results 
leads  to  the  conclusion  that  any  error  in  the  Dynamitron  energy  calibration  is  not 
statistically  significant. 
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ARL  Experiment. 


The  measured  activity  from  the  ARL  experiment  was  nearly  an  order  of 
magnitude  less  than  the  calculated  activity,  with  a  difference  in  excess  of  3  standard 
deviations: 

Theoretical:  Aq  =  42.92  ±  11.19  Bq 

Experimental:  A0  =  4.55  ±  0.04  Bq 

The  discrepancy  between  the  theoretical  and  measured  y-ray  activities  suggest  that 
either  the  current  or  the  energy  are  in  error  on  the  Varian  L200A  accelerator.  The 
possibility  that  the  energy  is  less  than  1.3  MeV  is  more  probable  than  the  current 
begin  in  error  by  a  factor  of  9.5. 

Based  on  the  results  of  this  simulation,  it  is  likely  that  the  estimate  of  1.3  MV 
for  the  Varian  operating  voltage  is  in  error,  and  the  actual  value  is  closer  to  the 
nominal  setting  of  1.0  MV.  Because  the  intermediate  state  at  1078  keV  has  an 
integrated  cross  section  16  times  larger  than  the  combined  cross  sections  at 
934  and  941  keV,  the  yield  of  115mIn  is  largely  due  to  the  1078  keV  intermediate 
state.  At  endpoint  energies  below  1078  keV,  irradiation  of  the  indium  sample  will 
produce  a  substantially  lower  yield  of  the  115mIn  isomer.  Thus,  a  likely  lower  bound 
on  the  operating  voltage  of  the  Varian  L200A  is  ~  1.1  MV. 

1.6  Conclusion 

The  results  of  the  simulation  and  the  WSU/AF1T  experiment  show  that  the 
theoretical  and  measured  y-ray  activities  due  to  the  photoexcitation  of  115mIn  agree 
to  within  a  standard  deviation,  so  there  is  no  statistically  significant  error  in  the 
Dynamitron  energy  calibration.  In  contrast,  the  significant  disagreement  between 
theoretical  and  measured  activities  resulting  from  the  ARL  experiment  is  evidence 


94 


that  suggests  the  Varian  L200A  energy  was  less  than  1.3  MeV.  At  a  nominal 
operating  voltage  of  1.0  MV,  the  results  of  the  GEANT4  simulation  suggest  that 
the  actual  operating  voltage  Vop  is  in  the  range  1.1  <  Vop  <1.3  MV.  It  is  suggested 
that  the  Compton-Hall  measurements  described  in  [50]  be  repeated  to  obtain  a 
more  accurate  estimate  of  the  actual  operating  voltage  for  the  Varian  LINAC. 

This  experiment  demonstrated  that  the  Dynamitron  at  WSU  is  capable  of 
achieving  photoexcitation  of  nuclear  isomers.  However,  its  limited  current  of  less 
than  1.5  /i A  at  operating  voltages  above  1  MV  represents  serious  problem  for  future 
experiments  involving  nuclear  isomers.  To  effectively  perform  nuclear  structure 
investigations  using  (7, 7')  reactions,  a  radiation  source  should  be  able  to  produce 
currents  of  >  10  /iA  up  to  at  least  1.2  MeV,  and  lower  currents  up  to  1.5  MeV.  An 
investment  in  radio  frequency  (RF)  tuning  the  Dynamitron  or  improving  the 
vacuum  system  may  help  achieve  these  gains.  In-depth  studies  of  nuclear  isomers  at 
WSU  would  also  greatly  benefit  from  the  installation  of  a  HPGe  cave-type  detector 
near  the  Dynamitron.  This  would  reduce  the  transit  time  between  irradiation  and 
counting,  which  is  essential  for  the  study  of  isomers  with  short  half  lives. 
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